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Conservation of energy implies that we cannot freely create superpositions between energy eigen- 
states, while access to coherence enables such processes. A simple example of the latter is when an 
atom is put into a superposition of two energy levels, using coherent radiation in form of a laser. The 
question is to what extent the coherence is degraded when utilized. Here we show that coherence 
can be turned into a catalytic resource, meaning that it can be used indefinitely to induce coherent 
operations without ever being degraded. As an application we consider quantum thermodynamics, 
where we demonstrate that access to sufficient coherence enables the release of all latent expected 
work content in a system. Since we can make it a catalyst, no coherence is spent to achieve this work 
extraction. In a wider perspective, the fact that coherence can be made catalytic provides a step 
towards a characterization of the relation between energy and coherence as resources in quantum 
statistical mechanics. 
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Introduction.- Coherence is a fundamental resource 
that enables us to implement coherent operations on 
quantum systems; a canonical example being the use 
of a laser to put an atom in superposition between two 
energy- levels, or radio-pulses for nuclear spins. Due to 
conservation of energy, we cannot directly turn a quan- 
tum system with a definite energy into a superposition of 
different energies. However, with access to a resource in 
form of a system with a high degree of coherence we can 
overcome this limitation to an arbitrarily good approxi- 
mation. The question is to what extent this 'coherence 
resource' gets degraded when utilized. 

Here we show that coherence can be turned into a cat- 
alyst, meaning that it can be used repeatedly without 
ever being degraded. We consider two models with this 
catalytic property. The first of these models is convenient 
to analyze, but is somewhat unphysical in that it has no 
ground state. The second model amends this problem, at 
the price of a protocol that regularly injects energy into 
the system. To touch upon the question of the existence 
of these effects in more realistic systems, we furthermore 
make a brief numerical analysis of the Jaynes-Cummings 
model [Hig. 

As an application we turn to the question of 'work ex- 
traction', i.e. how much work that can be extracted from 
a system in a non-equilibrium state with respect to a heat 
bath of given temperature. The work extraction problem, 
and the closely related concepts of information erasure 
and Maxwell's demon, have a long history (see e.g., [3]- 
|5]) with a recent revived interests, e.g., in the contexts 
of resource theories [5HTT] and single-shot statistical me- 
chanics j,12 18J. Recently, it was demonstrated [19] that 
for systems with non-classical states, meaning superpo- 
sitions between different energies, not all the expected 
work content can be extracted, if we are restricted to op- 
erate on the systems individually, and have no access to 
additional coherence resources. Here we show that access 
to sufficient coherence does enable the release of all the 
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latent expected work content of a system, even for indi- 
vidual operations. Since this analysis is performed within 
a model where coherence is catalytic, this implies that the 
additional release of expected work is not payed for by 
a loss of coherence. Finally, we make some observations 
in relation to single-shot work-extraction P^HT5] . which 
raise questions concerning the relation between coherence 
and noise-free energy. In relation to this thermodynamic 
application, one may note |20) . which used a coherent 
extraction device as a negentropy source, to demonstrate 
a transient efficiency that exceeds the standard Carnot 
bound for work extraction against a hot and a cold heat 
bath. 

Coherence as a resource. - Imagine a two-level system S 
with energy eigenstates |V'o) and IV'i); with corresponding 
energy eigenvalues hg and hi with hi > ho- If we wish to 
perform non-trivial operations on S, we need, due to con- 
servation of energy, another system where we can with- 
draw and deposit energy. We refer to this as the 'energy 
reservoir', E, with Hamiltonian He. If the spectra of Hs 
and He are suitably matched, i.e., if there are resonances, 
then the combined Hamiltonian Hs + He have degenera- 
cies. (For the sake of notational simplicity we avoid to 
write the more correct expression Hs ® 1e + ^s ^ He-) 
This means that the ground state \ipo) can be excited 
into IV'i) if we simultaneously de-excite the energy reser- 
voir. But what if we wish to transform the the ground 
state \tpn) into a superposition (jV'o) + |^i))/'\/2? With 
a bit of thought one can realize that it is impossible to 
do this with an energy preserving unitary operation if 
both the system and the reservoir initially are in an en- 
ergy eigenstate. (This observation can be understood in 
the much wider context of 'reference frames' and sym- 
metry preserving operations pTH^ .) At first sight this 
might seem a bit odd, since putting systems into super- 
positions between energy eigenstates is a common prac- 
tice in many labs. As is well known, one way to resolve 
this apparent paradox is to realize that these operations 
usually are induced by laser or radio-fields that can be 
modeled as coherent states, which are broad superposi- 
tions of the energy eigenstates (the number states) of a 
mode of the electromagnetic field [MH^B] (although this 



can be debated, see e.g. [27H32| ). Since coherence thus is 
a resource that allows us to perform non-trivial quantum 
operations, it appears a relevant question to ask to what 
extent it gets degraded when we use it. In the following 
we shall consider a model where it can be shown that the 
coherence does not degrade at all. In other words, it is a 
catalyst. 

The doubly-infinite energy-ladder.- The model we use 
has previously been considered in the context of quantum 
thermodynamics, for the purpose of work extraction [19]. 
The Hamiltonian of the reservoir in this model has the 
form 
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where s > is the energy spacing in this doubly infinite 
ladder, and {|i)}jez is an orthonormal basis, where Z 
denotes the set of (positive and negative) integers. Re- 

(s) 

garded as a Hamiltonian, _ff jj ' is slightly odd, in that it 
does not have any ground state. We shall shortly remedy 
this problem, but due to several convenient properties of 
this model (see Appendix A 1 1 we shall use it for the ini- 
tial analysis. We let S* be a two-level system with hi = s 
and ho = 0, and thus S and E are in resonance. (See Ap- 
pendix IA] for the case of A^-level systems.) Irrespective 
of which state the energy reservoir is in, it can always 
compensate for a transition on S. Next, we define the 
unitary operator A = J2jez, \J + l)(j|- ^s seen, this op- 
erator translates the state 'rigidly' along the energy lad- 
der. With the help of A we can formulate the following 
family of unitary operators on "Hs "SJ He 



A" 



ViU)= J2 Hn){^n\U\i>n'){^n'\- 
n,n''— 0,1 



(2) 



where U is an arbitrary unitary operator on Hs- As one 
can see, the mapping V does in some sense create an 
infinite number of copies of U on the combined space 
^S (X" "He- By construction, all V{U) commute with 
Hs -\- He, i.e., they are energy conserving (see Appendix 



Al). Furthermore they commute with all A°, and thus 
act uniformly over the energy ladder (see |19] for dis- 
cussions on this). The family of unitary operators V(U) 
serves as the set of 'allowed operations' in our model. In 
what follows, we shall investigate what kind of operations 
that this model allows us to perform on system S, and 
how this depends on the coherence of the state of the 
reservoir E. Finally, we show that this set of operations 
does not change, irrespective of how many times we use 
the reservoir. In this sense, the coherence of the reservoir 
is a catalytic resource in this model. 



FIG. 1: Generalized coherent states. As a simplified 
versions of coherent states, one can take uniform superposi- 
tions over a range of energy levels \rjL,io) = X]!=o \h) + l)/VT^- 
When interacting with a two-level system, the whole state gets 
shifted up or down along the ladder, depending on whether 
the two-level system absorbs or donates energy. For large 
L, the shifted and unshifted states have an overlap close to 
1, i.e., they are almost identical. This makes it possible to 
approximately create superpositions in the two-level system, 
and more generally to approximately perform unitary opera- 
tions. In the limit of large L, these implementations can be 
made perfect. 



Having defined our model, the question is what kind 
of operations it allows us to perform on S. Given a state 
a on the reservoir, we can implement channels on S via 



<^a,u{p) ^ TrE[V{U)p ® aV{U) 



(3) 



Given a state a we let C{a) denote the set of channels 
$(j_(7 that can be implemented via arbitrary unitary op- 
erators U. 

Next, we shall demonstrate that if a is sufficiently co- 
herent, then we can, to a good approximation, implement 
arbitrary unitary operations on S (see Fig. Tl). For this 
purpose, we consider states a — \riL.ia){VL.io\^^ the form 
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i.e., uniform superpositions over a collection of consecu- 
tive energy eigenstates. One can show (Appendix A 3 ) 
that in the limit of large L, the channel ^\til i ){vl i \,u 
converges to the unitary operation p i— > UpU^ . In other 
words, the higher degree of coherence (i.e., the larger L) 
the better we can implement arbitrary unitary operations 
on S. Hence, our model is powerful enough to achieve all 
unitary operations on the system, given a sufficient re- 
source of coherence in the reservoir. One may note that 
this is analogous to a similar constructions for coherent 
states of a bosonic mode interacting with an atom (see, 
e.g., chapter 7 in [33]). 



Catalytic coherence in the doubly-infinite ladder.- Let 
us now take a closer look on what aspects of the state a 
it is that determine the set of channels C{a) that can be 
reached. If one inserts Eq. (pi) into (pi), it turns out that 
^a-,u depends on a only via expectation values of the form 
Tr(A"cr) for a e Z. (See equation ( |A10[ ) in Appendix 



A2|) Hence C{cr) only depends on the set of numbers 
Tr(A"cr) (which does not determine a uniquely). 

To investigate the catalytic aspects of the coherence we 
need to determine how the state of the energy reservoir 
changes when we use it. For that purpose we define the 
corresponding channel on E 

Ap.u{<j)^TYs[V{U)p<E,aV{U)^. (5) 

By using Eq. (pi) one can confirm that 



Tr[A'^Ap,a(a)]-Tr(AV), 



(6) 



for all cr, p, U, and a. (See Lemmaplin Appendix |A4[ ) In 
other words, the expectation value of all powers of A are 
invariant under the action of these operations. Above, 
we noted that the only aspects of a that determine what 
channels can be implemented, are precisely these pow- 
ers (A°) = Tr{A°'a). We can thus conclude the set of 
channels that can be implemented the second time we 
use the reservoir is not affected by the first use, i.e., 
C{Ap^ij{c7)) = C{(j). ( For the full statement, see Propo- 
sition p] in Appendix A4 ) In other words, we do not 
degrade the coherence resource in the reservoir by using 
it. In this sense, coherence is catalytic in this model. 
One can also prove a stronger catalytic property, which 
does not assume that the reservoir and the subsequent 



systems initially are uncorrelated (see Appendix A 5 ) 

It should be emphasized that although (A") are invari- 
ants, the underlying state does change. Typically it ap- 
pears to get 'smeared out' over large r and larger stretches 
of the energy ladder. (See Appendix jA 6| for an example.) 
The half-infinite ladder (harmonic oscillator).- Since 
the state of the reservoir gets broader for each use, one 
could worry that the catalytic properties of the coherence 
in the above model is an anomaly, related to the fact that 
the model lacks a ground state; for a system with ground 
state the broadening distribution would eventually 'hit 
the bottom'. This is indeed true, but here we show that 
by introducing a very simple protocol, the set of channels 
that the reservoir can implement, can be be maintained 
indefinitely, also in a model that has a proper ground 
state. To this end, we 'cut away' the lower half of the 
doubly-infinite ladder and thus obtain (the spectrum of) 
the harmonic oscillator 
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We define a new class of unitary operations on S and E 
as 



-t-CX) 



with Vi as in equation (pi). By direct comparison one 
can see that V{U) and V+{U) act in an identical way on 
all states with at least one quanta, i.e, (?|cr|^) — for 
I = 0. (For multi-level systems this 'border zone' would 
be larger.) 

Let us now consider the dynamics in this model. For 
this purpose, assume that the state a have a projection 
only on a finite range of energy levels. Since the two- 
level system S can absorb or donate at most one energy 
quanta, it follows that the state on the energy reservoir 
can extend the support on the energy eigenstates with at 
most one energy level up as well as down in energy when 
the systems interact. Imagine now that reservoir initially 
is in a state with at least k quanta (i.e., (j|t|j) = for 
all J < k). For k uses of the reservoir on a sequence 
of two-level systems, the half-infinite ladder would thus 
be indistinguishable from the doubly-infinite ladder, and 
thus the set of channels that the reservoir can imple- 
ment remains constant for these k steps, i.e., a kind of 
'quasi-catalytic' property. By translating the initial state 
further and further up in energy (without increasing the 
'width' of the superposition), we can make the 'lifetime' 
of this quasi-catalytic behaviour arbitrarily long. (See 
Proposition H in Appendix |B 2| ) 

Protocol for a catalytic half-infinite ladder.- Next we 
shall see that we can achieve better than a mere quasi- 
catalytic property for the half-infinite ladder model. By a 
simple protocol we can maintain the coherence properties 
of the reservoir indefinitely (see Fig. p] and Appendix B 3 



for more details). To this end, we assume an initial state 
CTin such that {l\ain\l) = for / = 0,1, i.e., there are at 
least two quanta in this state. We let E interact with the 
first system Si with respect to some arbitrary choice of 
unitary y_^(C/i). As we know from the above reasoning, 
the effect is identical to V{Ui). We also know that the 
new state a on the energy reservoir is such that {1\(t\1) = 
for / = 0. Consider an ancillary two-level system A, 
with the the two eigenstates |ai), |ao), corresponding to 
energy the energies s and 0, respectively. We assume that 
A initially is in state |ai). By applying the operation 
V+{Ua) for Ua '■— \ao){ai\ -\- |ai)(ao|, we find that the 
state a is translated up one rung along the ladder, to a 
new state (Tout- Hence, the system is again in a state with 
at least two quanta. Since these operations all have been 
performed on states safely away from the ground state, 
we can conclude that C(tTout) = C(crin). By iterating this 
procedure we can thus conclude that the set of channels 
C that can be reach by this reservoir, can be maintained 
for arbitrarily many transformations. (See Proposition p^ 
in Appendix B 3 for more details.) 
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Decay of coherence in the Jaynes-Cummings model. - 
Although for many theoretical purposes it is enough to 
have demonstrated the possibility to make coherence into 
a catalytic resource, it is nevertheless relevant to ask to 
what extent these phenomena are noticeable in more gen- 
eral systems, especially if one would consider experimen- 
tal investigations of these questions. For this purpose, 
we consider the Jaynes-Cummings (JC) model [Il[l|, and 
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FIG. 2: Regenerative catalytic cycles. When the energy 
reservoir interacts with a two-level system (left part of the 
figure) the projection of its state onto the number states can 
increase with at most one level up and down in energy. As 
long as the projection onto the ground state of the reservoir 
remains zero, the set of channels that the reservoir can in- 
duce on the system stays intact from one interaction to the 
next. By using another two- level system in a pure excited 
state (right part of the figure) to inject energy into the reser- 
voir, the state can be translated 'rigidly' up along the energy 
ladder by one step. By alternating every use of the reservoir 
with such a pumping, the state of the reservoir can be kept 
away from the ground state, thus maintaining the coherence 
properties indefinitely. Note that the state of the energy reser- 
voir does change from one cycle to the next, e.g., the range of 
number states onto which it projects may become broader for 
each step. However, the relevant aspects of the state, which 
determines its capacity to induce channels, remains constant. 



provide numerical evidence that, akin to the half- infinite 
ladder model, the life-time (in the 'softer' sense of de- 
cay rates) of the coherence of a fixed- width superposition 
can be made arbitrarily long by increasing its average en- 
ergy (see Fig. [3] and Appendix O) . The JC Hamiltonian, 
H]G = g<^+ ^ a + gcf- (8) a^ , is a common model for in- 
teractions between a single electromagnetic mode and a 
two-level system. Here a, a^ are the standard annihilation 
and creation operators [a,a^] = l^, and u^ — |^i)(^o|, 
(J- = |7/'o)('0i|- As described in figure^ we measure, 
via the fidelity, the capacity of the reservoir to put a se- 
quence of two-level systems into superpositions of their 
energy levels. As seen, the decay of the fidelity appears 
to become slower the further up in energy the initial state 
is. (For more details, see Appendix C 2 ) 

Application: Coherence in expected work extraction. - 
In the work extraction setup we are given a system in a 
state characterized by a density operator p and a Hamil- 
tonian Hs- We are furthermore given access to a heat 



FIG. 3: 'Lifetime' of coherence in the Jaynes- 

Cummings model. These graphs depicts the decay of the 
fidelity by which the superposition 10) — (|?/)o)— i|'(/'i))/V2 can 
be created from the ground state |V'o), in a collection of two- 
level systems that sequentially interact with one single reser- 
voir, according to the JC model for a fixed time-step. Each 
curve corresponds to a different initial state of the reservoir, 
and shows the square of the fidelity i^*™' = (<;/)|<E>(cri?')|<^), 
against the number of times k the reservoir has been used. 
Each initial state o-^ = \rim){'rim\ of the reservoir is is 
a uniform superposition \rim) ~ X^;=_24 Ko('^^) + I), where 
lo(m) = (4m -|- 1)^, for m = 4, 6, 8, ... , 28, 30, counted from 
below on the right hand side of the figure. These graphs sug- 
gest that the lifetime (counted in number of iterations) of the 
coherence can be made arbitrarily long merely by increasing 
the energy of the state, without changing the 'width' of the 
initial superposition. This behaviour echoes the properties of 
the half-infinite ladder model. The dotted line is the value 
of (1 + |(7?„|A|7?„)|)/2 = 0.99, which would be the fidelity 
reached in the doubly-infinite ladder-model for these initial 
states. 



bath of temperature T, and our task is to extract as much 
useful energy as possible by equilibrating S with respect 
to the heat bath. Here we consider the question of how 
much work that can be extracted in an average sense, 
i.e., we strive to optimize the expected yield. Standard 
results O I34H36] suggest that the expected work content 
of a system is characterized by 



Asta„dard'(p,i?s) = kTD{p\\G{Hs)) , 



(9) 



where D{p\\r]) — Tr(plnp) — Tr(pln77) is the relative von- 
Neumann entropy, k is Boltzmann's constant, G(Hs) — 
exp{-l3Hs)/Z{Hs), Z{Hs) = Trexp(-/?ff5), and /? = 
l/(fcT). (With a base 2 logarithm, log2, in the relative 
entropy, we would get an extra factor ln(2) in ([9]).) In 
models that do not include an explicit energy reservoir 
one can confirm (see Appendix O) that the optimal ex- 
pected work content indeed can be characterized by equa- 
tion ([9|. 

However, in light of the above considerations on coher- 
ence, such approaches appears to implicitly assume the 
access to some form of ideal coherence resources. As men- 



tioned earlier, it has been shown [T^ that without access 
to any coherence, the optimal expected work content is 
not characterized by ([9]), but rather by 

^diagonal (P,ifs) = kTD([p]Hs\\GiHs)) , (10) 



where [p]hs — X); PipPu where Pi are the projectors onto 
the eigenspaces of Hg. The operation [■]hs i^ ^^ exam- 
ple of a 'pinching', where we here remove all off-diagonal 
blocks with respect to energy eigenspaces of the Hamilto- 
nian. However, one can show (Appendix p| that access 
to coherence increase the amount of work that can be 
extracted. Furthermore, in the limit of very large de- 
gree of coherence (large L for the energy reservoir states 
hi.io)); one regains equation (p]). (See Appendix E3 



From these observations we can conclude that coherence 
sails up as an important resource alongside with expected 
work content in a quantum thermodynamic setting. The 
question is how they relate to each other; how much of 
one resource can be gained by spending the other? The 
fact that we here perform the work-extraction analysis 
entirely within the doubly-infinite ladder model, implies 
that we only use the coherence catalytically. In other 
words, we do not 'spend' any coherence at all. 

Coherence and single-shot work extraction.- We have 
here focused on expected energy gain in work extraction. 
However, these results nevertheless suggest a couple of 
observations in relation to single-shot work extraction. 
As observed in the single-shot setting [T^18j . the ex- 
pected work content may not always correspond to or- 
dered 'work-like' energy (see discussions in [Til). Al- 
though the coherence is preserved in the dynamics of 
the doubly-infinite ladder-model, the state nevertheless 
appears to gradually get more broadly distributed over 
the ladder (see Appendix A 6 for a simple example). Due 



to this one may suspect that the 'quality' of the energy 
decreases. This observation maybe gets even more pro- 
nounced for the half-infinite ladder model, where we in- 
ject energy in the form of pure excited states to maintain 
the coherence. Since such pure excited states correspond 
to very ordered energy, the spreading of the distribution 
again suggest a degradation of the quality of the energy. 
(One may note that in terms of expected energy, the 
injection is not lost, but merely adds to the extracted 
energy.) These observations raise the question whether 
there is a cost in terms of ordered energy (as opposed to 
expected energy) associated with catalytic coherence. If 
one would compare models that do and do not preserve 
coherence, would there in some sense be different costs of 
noise-free energy for performing operations? Such inves- 
tigations would require the development of a (preferably 
operational) measure of ordered energy that takes into 
account coherence in the energy carrier per se. 

Conclusions and outlook.- Here we have shown that 
there exists models where coherence is a catalytic re- 
source. These models consists of a specific choice of 
Hamiltonian for a reservoir that carries energy and co- 
herence. Furthermore, they specify rules on the allowed 
dynamics in the interactions between a target-system and 



the reservoir. Here we have focused on the set of chan- 
nels that the reservoir can induce on the system via the 
allowed interactions, and demonstrated that the more co- 
herence there is in the reservoir, the more coherent the in- 
duced operations become. In the 'doubly-infinite ladder 
model' we have proved that the set of induced channels 
is invariant under repeated applications of the reservoir, 
and in this operational sense, the coherence resource in 
the reservoir is thus catalytic. To remedy the problem 
that the spectrum of this model is unbounded from be- 
low, we introduce the 'half- infinite ladder'. For certain 
classes of initial states of the reservoir, and a procedure 
where energy is regularly injected into the reservoir, this 
model simulates the doubly-infinite ladder, and is thus 
also catalytic. We have furthermore made a brief nu- 
merical analysis of the Jaynes-Cummings model, which 
indicates a behaviour that is reminiscent of the more ide- 
alized models. As an application of these findings, we 
have analyzed work extraction on non-equilibrium states 
against a fixed-temperature heat bath, and shown that 
a catalytic use of coherence can release all the expected 
work content of a system by individual operations. 



As mentioned above, we have here analyzed the coher- 
ence in the reservoir in terms of the channels that the 
reservoir can induce. For the models we have considered, 
the 'relevant' aspects of the state of the reservoir are 
the expectation values (A°). This certainly depends on 
the 'allowed dynamics' we have declared in the models, 
and thus different models could potentially single out dif- 
ferent 'types' of coherence as 'relevant', and potentially 
invariant. How this should be understood in a wider per- 
spective requires further analysis. As maybe already has 
been indicated, the question of coherence is closely re- 
lated to the concept of reference frames and symmetry 
preserving operations pTM23j . The considerable machin- 
ery that has been developed in this context may favorably 
be applied to further analyze the question of catalytic co- 
herence, but also to provide alternative ways to analyze 
these questions. 



We have here argued that coherence is an important 
thermodynamic quantity. One may thus ask why we 
normally do not encounter coherence in discussions on 
thermodynamics. In the limit of many independent sub- 
systems in identical states, and identical non-interacting 
Hamiltonians, the global state has a low degree of 'off- 
diagonality' with respect to the joint energy eigenspaces 
[TUl [T^ [TU] and hence the difference between (|9| and 
( 10 1 becomes small. Since many realistic systems can 



be argued to be 'close' to this situation, we would thus 
normally not expect to see the effects of coherence. How- 
ever, as demonstrated above, this is no longer true when 
we step outside of this 'multi-copy' setting. To obtain 
a complete theory of quantum thermodynamics, this is 
arguably a necessary step. 
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Appendix A: Doubly infinite energy ladder 

As a first step to investigate the catalytic properties 
of coherence, we here consider a simple model for en- 
ergy conserving operations on a combination of a system 
and an energy reservoir. This model has previously been 
considered in \Wi. 



1. The doubly-infinite ladder model 

Our model consists of two parts: A 'system' S, mod- 
eled by a finite-dimensional complex Hilbert space H5, 
and an 'energy reservoir' E^ modeled by a separable 
Hilbert space He- We let {|j)}jgz denote a complete 
orthonormal basis of He- Here Z denote the set of (neg- 
ative as well as positive) integers. Given a real number 
s > we define the following unbounded Hermitian op- 
erator 



m 



^E^1^')01- 



(Al) 



We will in the following regard this operator as the 
Hamiltonian of the energy reservoir E, although, needless 
to say, it is somewhat unphysical, since it has a 'bottom- 
less' spectrum. However, we will remedy this problem is 
Section p] ( |19j mostly uses a continuum reservoir rather 
than a discrete model. Although see the remarks in Ap- 
pendix G of [H].) 

Given a finite-dimensional Hilbert space Hs and a real 
number s > we let Hs{'Hs) denote the set of Hermitian 
operators on %$ such that all its eigenvalues are integer 
multiples of s, i.e., all its eigenvalues can be written as 
hn — szn where z„ € Z. The reason for this construction 
is to 'match' all energy differences in Hg with the energy 

is) 

spacings in H^ . By this we make sure that for each pair 
of eigenstates in S there is a transition (in fact infinitely 
many) in E that can compensate for the change in energy. 

The eigenspaces of the the joint system Hs + H^' (we 



should strictly speaking write Hg ( 
a very simple structure. 



i£; + is«)ijj;''^)have 



Lemma 1. Let s > be given. Let Hs be a finite- 
dimensional Hilbert space with dimHs = N, Hs G 
Hs(Hs), with eigenvectors {ipn) o-nd corresponding eigen- 



values sZn (where z„ G ZJ. Then 



Hs + H^E^ ^s 



3& 

N 



(A2) 



pW:=^|V>„)(V„|0|j-^„)(j 



n=l 



is an eigenvalue decomposition, where P^^' is the projec- 
tor onto the eigenspace corresponding to the eigenvalue 
js- 

The whole purpose of the construction with the en- 
ergy reservoir is to keep track of energy. Hence, for the 
dynamics on the joint system SE we consider unitary 
operations that conserve the total energy, i.e, that com- 
mute with the joint Hamiltonian Hs-\-H^ . However, we 
furthermore demand that the dynamics should act 'uni- 
formly' over the energy levels in the reservoir. To express 
this a bit more precisely, we introduce a unitary operator 
A on He, which can be thought of as a 'rigid translation' 
along the doubly-infinite energy ladder. 



A:=E|A:-M)(fc|. 



kei 



One can see that 

A" A'' = A'' A" = A''+^ 
where we let A" = Ik. 



A' 



at 



A- 



(A3) 



a, 6 G Z, (A4) 



Definition 1. Let s > 0, and let Hs be a finite- 
dimensional Hilbert space and let Hs G Hs{Hs)- We 
use the following notation 

• L{Hs) is the set of linear operators on Hs, and 
U{Hs) the set of unitary operators on Hs- 

• \Jf{Hs+HE) is the set of bounded linear operators 
X on Hs (Xi He such that [X, Hs + H^ ] — 0, and 
[X, is (g) A"] = for all a G Z. 

• VT{Hs + H^ ) is the set of unitary operators U 

on Hs (X" He such that [U, Hs + H^e ] ~ ^' '^^'^ 
[U, is (g) A°] = for all a G Z. 

• Given a Hilbert space H we let S(H) denote the set 
of positive semi- definite trace class operators with 
trace 1 onH, which we refer to as density operators 
or 'states'. We let S+(H) denote the set of density 
operators with full support. 



The set \JT{Hs + H^ ) is precisely the set of evolu- 
tion operators we allow in our model; they conserve en- 
ergy, and commute with all energy energy translations 
A". The following lemma shows that we can describe 
the set UT(i/5 4- H^ ) in a very convenient manner. 

Lemma 2. Let Us be finite- dimensional with N :~ 
dim'Hg, and let Hg S Hg{H.g). Let \'4>n) be eigenvec- 
tors of Hs with corresponding eigenvalues sz„ . Define 
the mapping 

N 
V{Q):^ Y, |V'n)(^n|g|V'n')(V'n'|<^A^"'-^", (A5) 



for all Q in L{Hs). The mapping V is a bijection between 

L{T-Ls) andhT{Hs+H^ ), which preserves the structures 
of these operator spaces, in the sense that 



V{aA + I3B) = aV{A) + I3V{B), 
V{AB) = V{A)V{B), 
V{A^) = V{A)\ 



for all A, B G L{Hs) and all a, /3 e C. As a consequence, 
V is a bijection between U{'Hs) and \]T{Hs -\- H^'). 



Equation ( A5 ) provides us with a very convenient han- 
dle on VT{Hs-{~H^'), since we can reach each element of 

\]T{Hs -h H^E^) uniquely by V{U) for some U e C/(Hs)- 
Due to this reason we will in this investigation hardly 
ever refer explicitly to the set 1]T{Hs-\-He ), but rather 
describe it indirectly via V. 

It is worth noting that one can rewrite the mapping V 
as 



V{Q) ^J2^,{Q), 



N 

VjiQ) - Yl IV^«)(V'"l<3IV'n')(V'n'| <» \j - Zn){3 - Z„'|, 
n,n' — l 

(A7) 

where V,(Q) = P'^^WjiQ) = Vj{Q)P^^\ Hence, what V 
essentially does is to embed an infinite number of 'copies' 
of Q into the space of operators on T-Ls ® T~Le- Moreover, 
if [/ e U{Hs) then Vj{U)Vj{UY = Vj{U)^Vj{U) = P'^''. 

Proof of Lemma^ Since V{Q) merely is a finite sum of 
bounded operators, it follows that V maps LiT-Lg) into 
the set of bounded linear operators on T-Ls ® 'He- Since 
all A" commutes, it follows that [V{Q), Is ® /S.'^] =0 
for all a € Z. One can also verify (e.g. via showing 

[H^E^^''] = saA") that [V{Q),Hs + H*^^] = 0. Hence 
V{L{ns)) C hT{Hs + hI^^). It furthermore only re- 



quires a verifications to see that the conditions in ( A6 ) 
hold. It thus only remains to show that V is surjective 



[i.e., V{L{Hs)) = LT(ffs + h'j^^)] and injective (i.e., 
two different elements in L{Hs) cannot get mapped to 
the same element). The injectivity can be obtained if 
we apply ( [AS] ) to V{Qi) = V{Q2) and note that A"" is 
never the zero operator. For the surjectivity we assume 



r(s) 



(s)i 



Y e LT(iJs + i?ijO> and note that [Y,Hs + H)^'] = im- 
plies that Y = Ejez-P'-'-'^-P^^^ with pf-?' as in Lemma 
[T] Hence, 



Y = 



EE 

j 71, 7l' 



{{^m\{j - Zr.\)Y{\i;^,)\j - z^,)) 



(A8) 



|V'«)(V^n'|® \j - Z„){j 



Using \j-zn,) = A^l-z^,), and [ls®A^']Y[ls(g>A^] = Y 
(since [Y, I5 (g) A-'] = 0) we can conclude that 

X ■■= J2(i^'n\{-Zn\)Y{\,Pn')\ " ;2„')) |^n) (^n' | , (A9) 

is such that Y = V{X), which proves the Lemma. Com- 
bining the bijectivity of V with (A6) yields that V also 

n 



(A6) i^ ^ bijection between U{l-Ls) and i]li{Hs + H 



(s)n 



2. Induced channels 

Given the unitary operators V{U) on the joint system 
SE we can induce channels on S and E. 

For a fixed state a on the energy reservoir, each choice 



of unitary U ox\ S induces a channel ^^ tj on system S 



via 



-^S,H, 



K:u^P)-^'^'^EV{U)p®aV{U)^ 

{tpn\U\'lpn'){'^n'\p\i^m'){'>}^m'\U'^\^ri 



N 



n,n',m,m' — 1 



|V'„)(V'm|Tr(A^"' 



-Zj^ — z I -\-z„ 



a). 



(AlO) 

We denote the set of channels that can be reached on S, 
by using cr as a resource, by 

C^^"{a):=[^l^^:UeU{Hs)}- (All) 

As seen, both '^^'jj and C^'^ depend on the choice of 
Hamiltonian Hs, both via the eigenstates \^n), and the 
eigenvalues sZn- However, we will in the following of- 
ten drop one or both of the superscripts when they are 
not needed, i.e., <&^'[/ , ^^c/i ^"^^ ^(y,u denote the same 
object. Similarly for C^'" , C^ , and C. 

We can also express the effect of V{U) on the energy 
reservoir via the channel 

Ap^ui^j) ■.=TTsViU)p®aViUy 



N 

E 



n.n' ,m' — l 

A^"'-^"(tA^'"'^ 



(A12) 



3. Approximation of energy-mixing unitary 
operations 

At first sight it may not be clear what kind of oper- 
ations one can reach on S. Here we show that in the 
hmit of large coherence, all unitary operators can be im- 
plemented on S. As mentioned in the main text, this 
model can be viewed as a simplified version of transitions 
in atoms induced by coherent electromagnetic fields (see 
e.g. Chapter 7 in [33]). In both cases the key-point is 
that the state of the energy reservoir is such that adding 
or removing a small amount of energy does not change it 
much. 

Define the states 



^N 



1 



\VL,lo 



L-1 
^ 1=0 



l + lo), 



(A13) 



i.e., this is a uniform superposition of L > 1 sequential 
energy eigenstates, starting at Iq. One can verify that 



(r7i,,JA"|77i^,„)=max(0,l-^) 



(AM) 



In other words, the shifted A"|?7L.io) and the unshifted 
state \ijL.io) ^re almost indistinguishable if L 3> a. 

In the following, we shall sec that the difference be- 
tween the induced channel ^a-,u a-nd the unitary channel 



Uu{p):^UpU\ 



(A15) 



becomes very small when L ^ a. Hence, in the limit of 
large L, we can implement arbitrary unitary operations 
onS*. 

To this end we do, on the set of super-operators 
r : L{'Hs) ~^ L{'Hs)j consider the operator norm in- 
duced by the trace norm ||r||tr := supuxn^^i ||r(X)||i. 
By trivially extending the super-operator to a sufficiently 
large ancillary Hilbert space Ha (it is sufficient with 
dim Ha — dim'Hs ,37|) one can define the diamond norm 
[371135] ||r||o := ||r®/A||tr, where Ia is the identity map 
on LCHa)- 

In the following, we let C"^" denote the set of n x n 
matrices with complex entries. Furthermore, for A, B Cz 
(£nxn ^g Yet A * B denote the Hadamard product, i.e. 
the component-wise product (A * B)jk '■= AjkBjk- We 
furthermore let || • || denote the spectral norm \\Q\\ :— 
supy^ll^i IIQV'II- Fo^ an arbitrary Q G £Nxn ^g jg^ j^^ 
the following {si{Qj)fLi denote the singular values of Q 
(including zero values) ordered non-increasingly si{Q) > 
S2{Q) > ■ ■ ■ > sn{Q) 



Lemma 3 (Theorem 5.5.4 on p. 334 in 
C^'^^ then 



J2s,iA,B)<J2^iiA)si{B), 



Let A,B e 



(A16) 



for k — 1, 



1=1 



,N. 



1=1 



Since llQIli = Eili^iiQ) and |1Q|1 = si(Q) it follows 
directly from the above lemma that 

P*i3||i<||A||||i3||i. (A17) 

Lemma 4 (Eq. (3.7.2) on p. 223 in [35]). Let C G C"^", 

then 



||C||< /maxJ^IQ,! /maxJ^IQ,! 



(A18) 



Proposition 1. Let s > 0, Hs with dim'Hs — N, and 
Hg e Hg{'Hs): with eigenvalues szn and corresponding 
orthonormal eigenvectors \4'n) ■ Let ^^^,11 be as defined in 
equation {AlO), and let |?7L,io) be as in equation {A13), 
for Iq € TL and i G N. Let z„ 
min„z„. // 



:= max,, z„ ana z- 



L > 2{z 
then for very U G U{%s), 






niin ■ 

(A19) 



and 



ll^y-'J'h.,,o>('7..„l,C/||o<2iV 



3 '-■max 



L 



(A20) 



(A21) 



Proof. Let T-La be a Hilbert space with orthonormal basis 
{|a„)}*£]^. (The choice of M will in the end give us the 
induced norm [M — 1) or the diamond norm (Af — N).) 
Let X G L{Hs ® Ua) be such that |lX|li = 1. For the 
sake of compactness we will use the notation \4'ki0.n) '■— 
\ipk)\an)- 

Due to L > 2(2:niax ~ 2niin) it foUowS that \zk — Zk' 



zi + 2i'| < L. Combined with equations (AlO) and (A14) 
this yields 

[U ® iA]X[U^ ® Ia] ~ [$f,,,,„)(„,,,„|,t/ ® Ia]{X) 

= J2 Y.(^k'\^\^k){->Pk,an\X\iJ;i,an'){lbi\U^\^i') 



k.k'Ll' n,n' 



|V'fc',a„)(V'z',a„ 



\Zk - Zk' - zi + Zl,\ 

L 



(A22) 



By the triangle inequality for the trace norm 

[U ® iA]X[U^ ® Ia] - [<I'f,,,„)(„,,,j,a ® IaKX) 

n^n' k,l 

(A23) 

where '*' denotes the Hadamard product, and where 
the matrices yl(fe'''«^"') := [^i^,;"'"'^]fc',/' and SC^'^') := 
[L^k''i']k',i' are defined as 



4^;;"^"'):=(V;fe,|[/|Vfc)(V;fe,a„|X|V^,a„,)(7/.^|f/t|V;^) 

(A24) 



„(fc,0 ,_ \zk - Zk' - zi + zi,\ 
k',i' — ^ 
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By Lemma ji) ||B('=^')|1 < N{z,^,^ - z^i^)/L. One can 
confirm p(fc.'."."')||^ = KV'fe, anjXlV'i, a„'>l- By combin- 
ing these observations with equations (A17) and (A23) 
we can conclude that 



[U ® 1a]X[U^ ® Ia] - K,.,^n,,,^iu ® Ia]{X) 



<N 



■^ V-2^max -^^niin } 



L 



^^\{'>i^k,an\X\'ipi,an, 



<MN 



n.n' k,l 
2 ■^V'^inax -^min j 

L ' 



(A25) 



where the last inequality follows by Cauchy-Schwarz. If 
we choose the ancillary Hilbert space to be trivial, i.e., 
dimHyi = M = 1, then we obtain the bound on the 
norm || • ||tr- It is sufficient to choose the ancillary Hilbert 
space to have the same dimension as "Hs to obtain the 
diamond norm [37]. In this case M = N and we obtain 
the corresponding bound. D 



4. A catalytic property 

The previous section shows that the set of channels 
C{a) which we can reach on system S depends on the co- 
herence in the state a of the energy reservoir E. In other 
words, the coherence in the state of £^ is a resource for 
the implementation of channels on S. Here we show that 
this resource is catalytic, i.e., if we use E to implement 
a channel on one system 5*1, then this does not diminish 
the capacity of E to subsequently implement channels on 
another system 82- It is worth noting that this catalytic 
property holds for all states a; it is not limited to states 
with a large degree of coherence. 

Consider the expression for the channel <I>^ ^ in equa- 



tion (|A10|). As one can see, the only way in which ^fjj 



depends on cr, is via the expectation value of powers 
of A, i.e., on the sequence of numbers (Tr(A°cr))^g2. 
(Strictly speaking, it only depends on a finite number of 
them.) One should note that this set of numbers does not 
uniquely define a. If one think of a represented as a ma- 
trix in the energy eigenbasis, i.e., in terms of the matrix 
[{jW\j')]j,j', one can see that Tr(A°CT) corresponds to the 
trace, Tr{A^a) = J^iUWlJ + 1) ^^ the sum of the ele- 
ments in the first upper diagonal, Tr(A^cr) to the second, 
etc. The following lemma shows that these expectation 
values are invariant under the action of the channel Ap^o- 
on the energy reservoir. 

Lemma 5. Let s > 0, Us finite-dimensional, and Hs G 
HsiHs). Then 



Tr(A'^Ap,c/M)=Tr(AV) 



(A26) 



for all a G S{He), U G UiUs), P e SiUs), and a G Z, 
where ^p,u is as defined in equation (A12). 



Proof. Let szk and jV'fc) be the eigenvalues and cor- 
responding eigenvectors of Hs, and put Uk' k '■= 

{il^k'\U\'il;k)- Then 

= J2 Uk',kU*^ri{^k\pHi)T^r{A''A'''-'>''aA''-'^'^) 

k,l,k' 

= ^ Uk',kU*k,^i{^k\p\^i)Tr{A-^'+^-'A^A^''-^^'a) 

k,l,k' 

= ^ Uk'.kU*k,A^k\pHi)Tr{A-+'^-^'a) 

k,l,k' 

-5]4,/(V'fc|p|V'/)Tr(A'^+^'=--V) 

k,l 

= ^(V'fc|p|^fc)Tr(AV) 

k 

=Tr(A''CT). 

(A27) 

n 

The following proposition tells us that if we use the 
reservoir twice, then the set of operations we can reach 
in the second application is independent of what we did 
in the first application. In other words, it is as if the first 
application had not happened. It is maybe also worth to 
point out that that the only assumptions we impose on 
the Hamiltonians of the subsequent subsystems are that 
they arc non-interacting, and satisfy the 'energy match- 
ing condition'. We do not assume that they are identical 
(or even operating on Hilbert spaces of the same dimen- 
sion). It is straightforward to generalize this proposition 
to more than two subsystems. 

Proposition 2 (Catalytic coherence). Let s > 0, and 
let T-Lsi and 'Hs2 ^^ finite- dimensional Hilbert spaces. Let 
Hs, G H.ins,) and Hg, G HsCHs,). Let Ui G C/(HsJ 
and pi G S{Hsi), ^2 G [/(Hsj) and p2 G ^('Hg^) . Le t 
Api,(7i be the channel on E as defined by equation (A12) 



for interaction with system Si, and let ^ '^ be the channel 
as defined by equation ^AlO ), for system 82- Then 

Sn.H 



^A^"uAa),uM) 



$; 



Equivalently, this can be reformulated as 
C^-^-KAp„t/,(a))=C^-^-^(a). 

Proof. Let IV'n ) and szn be orthonormal eigenvectors 
and corresponding eigenvalues of Hs,, and let |V'm ) 
and szm be orthonormal eigenvectors and correspond- 



(A28) 



(A29) 



ing eigenvalues of -ffgg. By equation (AlO) we can write 



$ 



S2,Hso 






k,k'.Ll' 



J2) (2) ^(2)_^^(2) 



l^(^))(4^'|Tr[A^r'-;;:'-r'+4-'Ap,^c/.(^)]. 



(A30) 
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By Lemma [5] we know that 



(2) (2) ^(2) , ^(2) 



^< +^<' A„, r;, (a)] 



(2) (2) (2) , (2) 



which proves (A28) and (A29l 



(A31) 



D 



In the above proposition we made no assumptions on 
the relation between the first and the second application 
of the reservoir. For the sake of clarity, we here also treat 
the special case that we do use the same Hamiltonian on 
both the first and second system. 

Corollary 1. Given the assumptions of Proposition^ 
hut with the additional assumptions "Hs^ ~ ^Si; H :— 
Hs2 — Hsi • *^ follows that 



C 



S2M 



(Ap,,t/.M)-C^-^(a). 



(A32) 



Strictly speaking the correct formulation of (A32) 
would be to say that the two sets of channels, 
C'^2'^(Apj_[/j (ct)) and C^^'^{a), are isomorphic (since the 
channels act on different spaces). 



Proof. Assume that Hs2 
with the isomorphism X 



and Hsi are isomorphic. 



,7V 



/,(2)\/./,(l) 



with z„ 



SZ' 



(1) 



\^i'') 



E„=i IV'n )(V'n I, where 



J2) 



and Hs2\->Jrn 



(2)\ 



SZ 



(2)l/,(2) 



M) 



\rn' 



Hsi^)- By equations (ABO) and (A31) we see that 
-S2.HS, 



(and thus '-ffs^ 



X^%:::M,uM)^ - K^wlxi^'P^^)- Hence, we 
can implement 'the same' channels on 5*2 as we can on 
Si, and thus, in this sense, equation (A32) holds. D 



5. A stronger catalytic property 

In the previous section we assumed that the system 
and the energy reservoir initially are uncorrelated (which 
indeed is necessary if we wish to the express the evolution 
in terms of channels). Here wc show that even if we 
initially have an arbitrary joint state on the reservoir and 
the systems, a catalytic property still holds. 

Lemma 6. Let s > 0, and let Hsi o,''^d 'Hs2 ^2 finite- 
dimensional Hilbert spaces. Let Hsi € HgCHsi) and 
Hs2 e Hs{Hs2), and define Hs := Hs^ ® Hs2 and 
Hs '■— I2 <8) Hsi + Hs2 ® Isi • I'he mapping V , as defined 
in Lemma [^ can he defined on the spaces T-Lsi ® ^E' 
'Hs2 fX" "HI;, and Us (8) H|;, with respect to the Hamilto- 
nians Hsi, Hs2, and Hs, respectively. For the sake of 
clarity we denote these with superscripts. Lt is the case 
that 



v"^{is2 



is J 



ls2®V"'HQi), ^QieLiUs,), 
V"'-{Q2)®is,, VQ2eL(HsJ. 



(A34) 



Furthermore 

[is2®V"'^{Qi)][V"'-{Q2)®is^] 

= [V"'- (Q2) ® isj[is2 ® ^"""^ (Qi)], 

for all Qi e L{Hsi) '^^d. Q2 € L{Hs2)- 

Proof. Let \4'n ) and szn be orthonormal eigenvectors 

(2) 

and corresponding eigenvalues of Hsi, and let \^jjm} 

(2) 
and szm be orthonormal eigenvectors and correspond- 
ing eigenvalues of i/g, . For a compact notation we de- 
fine |i/'m\'0n ) := |V'm^)|V'n )■ By usiug equation ( A5 ) 
we can write 

v"'iis2(®Qi) 

= E ('/'L'\^i')|(is.0Qi)IV'L'\V'iV) 



LS2 



v^^^^(Qi). 



An analogous reasoning holds for V^^{Q2 fX" Isi)- By 
Lemma [2] we know that V preserves operator multiplica- 
tion, and thus 

[is2®V''^^{Qi)][V"''^^^{Q2)®isA 

^v"'{is2(®QiWHQ2(^is,) 

^V"'{Q2®Qi) 

= V"^{Q2®ls,)V"'{is2®Qi) 



[V"'-{Q2)®is^][ls2®V"'^{Ql)]. 



u 



Suppose that we, apart from the energy reservoir E 
have two systems of interest ^i and S'2, and that the 
joint system S2S1E initially is in some arbitrary joint 
state rj. 

One can compare two different scenarios. In the first 
scenario, we first implement a unitary operation V^^^ on 
ESi followed by a unitary operation V^^^ on £^5*2 . In the 
second scenario, we only perform the second operation 
yHs2 on ES2. The following proposition tells us that 
the reduced density operator on ^2 is the same in both 
scenarios. 

Proposition 3 (Strong catalytic property). Let s > Q, 
and let Hsi and T^Sj be finite- dimensional Hilbert spaces. 
Let Hs^ € Hs{Hsi) and Hs^ G Hs{'Hs2), and define 
Us = "Hsi ® 'Hs2 and Hs = 12® Hs, -\- Hs2 ® Isi • Let 
Ui e U{ns,) and U2 e U{ns2), ihen 



Trs5, [V"'- {U2)V"-^ [UiW"^^ {Ui)^V"-^ (t/2)^ 
= TrB(F^^2([/2)Trs,(r;)y^-2([/2)t), 



(A35) 



(A33) for all r] e S{Hs, ® Hs2 ® ^-e)- 
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Strictly speaking we should rather write '[V-^^2 (JJ2) (g) 
isj[is2 <^ V"^i{Ui)y than 'V"^2{U2)V"^i{Uiy in the 
above proposition. 

Proof. By Lemma [6] and cyclic permutation under the 
(partial) trace Tr^Sj^ it follows that 



Tr^si 

Tr^si 



This provides an example of the spreading of the state 
mentioned several times in the main text. As a special 
case one can choose ct(°) := \vL,i„){VLdo\, for I'/jlAo) = 



=Tr 



ESi 



(V"'2 ([/2)F^si ([/^)^yffsi ([/^)t-5/Hs, (^2)t) 
(v"'^ {Ui)V"''2 {U2)vV"'^ {U2)^V"''^ (t/i)^) 
([is. ® V"'^ ([/i) V^«i (C/i)]F^-^ ((72)77!/^^^ ((72)^ 



(A36) 



According to Lemma [2] V preserves Hermitian con- 
jugate, operator multiplication, and the identity 
V^^i {Ui)'^V^^i (Ui) — isi- From this we can conclude 
equation ( |A35[ ). D 



X]/=o I' + Io)/VL. In this case we obtain F^ 



{k) 



1/{2L). 



opt 



= 1 - 



Appendix B: Half-infinite ladder (harmonic 
oscillator) 

As mentioned earlier, a problem with the doubly- 
infinite energy ladder is that it is somewhat unphysical, 
as it has no ground state energy. However, by 'cutting 
away' the lower half of the spectrum we get the spec- 
trum of the harmonic oscillator. Here we show that we 
can, with some modifications, reconstruct the catalytic 
properties of the doubly- infinite energy ladder in this har- 
monic oscillator model. 



6. An example 



The half-infinite ladder model 



Here we consider a simple illustration of the cat- 
alytic implementation of channels. Consider a collection 
of non-interacting two-level systems. We let \'ipo),\ipi) 
denote the eigenvectors with corresponding eigenvalues 
Zq = 0,zi = 1 of each system Hamiltonian. Suppose that 
all two-level systems initially are in the ground states 



\ipo), then the channels (AlO) and (A12| take the form 



*.,y(iV'o)(V'oi) ^\uoo\^m{^o\ + it/ion^i)(Vi 

+ C/rot/ooTr(Aa)|V'o)(V'i| 
+ C/ioC/o*oTr(Ata)|^i)(V'ol, 
^IVo>(^ol,t/(^) =\Uoo\^cr + |t/iopAVA, 



Here we use an harmonic oscillator as model (or any 
Hamiltonian that is iso-spectral to the harmonic oscilla- 
tor). 



+CX3 



H^--=sJ2j\j){j\. 
J=0 



(Bl) 



As before, we define an iV-level system with a Hamilto- 
nian Hs S Hs{'Hs), with eigenvalues /i„ = sz„. We also 
let 



min z„, Zmax := max z„. (B2) 

= 1....,N n=l,....N 



where [/„„' = {^n\U\^n')- 

Suppose we ideally would wish to put all the two- 
level systems into the state \(j)) = {\ipo) — »|^i))/'\/2, 
where we sequentially use one and the same reservoir, 
initially in state cr'"'. As a measure of the quality 
of this preparation we take the (square of the) fidelity 
Fif :— ((/)|$g.(o) u{\tPo){tPo\)\(I>), and optimize over all uni- 
tary U. An optimizing unitary is obtained by C/10 = 
l/\/2 and Uoo = -iexp[-iarg(Tr(Acr('^)))], and leads 
to the optimal value F^'^l = 5 + i|Tr(Ao-(o))|. with 
this optimizing unitary, the next state of the reservoir is 
(t(i) :== ifT(°) + iAV(°)A. As seen (and as we already 
know) Tr(Acr*^^'') — Tr(A(T*^°-'). Hence, in the second step 
we can use the same optimizing unitary, resulting in the 
same optimum F^pj = ^opt- Clearly, this process can be 
repeated indefinitely, with a constant sequence of fideli- 
ties F^pl = \ -h i|Tr(Acr(o))| for aU k. However, the state 
of the reservoir changes, and at the fcth step it is 



,(fe) _ 






A'V(°)A'. 



(A37) 



The projectors onto the eigenspaces of Rg + H^ can be 
written 

Pf-= Y. |^n)(V'n|®|^-^n)(;~^„|, (B3) 

)l=l,...,A':/>z„ 



for all I > 



(It is maybe slightly annoying that 



we let the index / start at Zmin rather than at zero, but 
we do this for the sake of compactness of formulas and 
coordination with the previous section.) Note that 



P« = 



Y, \'^n){'^n\<»\l-Zn}{l-Zr,\, VI > 



n=l....,N 

(B4) 
In other words, P^ = P^^^ for I > Zmax- (We will in 
this section often compare operators on the spaces HJ^ 
and He 7 where the former can be regarded as a subspace 
of the latter. Without further comments we do in these 
cases assume that the operators on "H^ are extended triv- 
ially as zero operators to the orthogonal complement of 
Hj^ in He, i-e., they act trivially on the negative half- 
ladder.) 



For every U € U{%s) we define the unitary operator 

Zmax-l 



XiX}^x}Xi^P^_l\ 



where 



N 



Vl{Q):= Y^ \i^n){i^n\Q\i^n'){'^n'\®\l-Z.a){l-Zr,>\, 
n,n' — 1 

(B6) 

is as in equation ( |A7| ). The unitarity of V^(\J) follows 
from the fact that Vi(f/)V^(f/)t = V^(C/)Vz(t/) = Pf. 
We will not pay any particular attention to the choice of 
the operators X/, as we will focus on states outside the 
supports of these operators. As seen, we have, by the 
very construction, made sure that y+(C/) in (A7) and 



y(f/) in (B5) act identically for states with support suf- 



ficiently far away from the ground state. The rest of this 
section is devoted to formalize the simple idea that if we 
start of the energy reservoir in a state that is sufficiently 
far from the ground state, we can maintain it so by suit- 
able arrangements, and by doing this we retrieve all the 
relevant properties of the doubly-infinite ladder model. 
Analogous to A we define the operator 



+00 



A+:=^b- + l)0-|. 



(B7) 



i=o 



Note that A_|_ is not unitary, but rather a partial isometry 
that takes H^ into the subspace Sp{|j)}j>i. (This is of 
course possible only because H^ is infinite-dimensional.) 



2. Sufficiently far from tlie ground state, tiie two 
models are equivalent 



Analogous to the channel ^ajj in Section A 2 we define 
*^,r/(p) := TrsK,(C/)p® ay+(C/)t. (B8) 

Similarly, we denote the set of channels that can be 
reached on 5, using cr as a resource, by 



C^^^^(a):={<l>+^:C/eC/(H5)}. 



(B9) 



The effect of V+(C/) on the energy reservoir is expressed 
via the channel 

A+c/(^):-Trsl^+(f/)p®ay+(C/)t. (BIO) 

Define the family of projectors 

^>rn:=^b')(j|, ™>0. (Bll) 

i>m 
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Lemma 7. Let s > 0, "Hs finite- dimensional, and 
Hs € Hs{'Hs) with maximal eigenvalue sz^ax cind min- 
imal eigenvalue sz„^in- Let a G S{'He) be such that 
P>z . -z ■ <^P>z -2, . — C; then 






(B12) 



for all U e UiUs). 



,(') 



Proof. By comparing with the definition of P\_ in (|B3|) 
one can see that 



j«ri 



0, 



nin ^ ^ "^ 



(B13) 



PX'[^S®P>z_______ 

Analogously 

Note that P|'^ = P(') and Vi(C/)P(') = VijU) for / > 



(B14) 



j(0 



Zjnax, and XiP\' = Xi. By comparing (A5), (AT), and 



( B5 ) , it follows that 



v+{u)[is®P> 



v{u)[is®P>z„ 



(B15) 
The statements of the lemma follows directly from this. 

a 

Lemma 8. Let s > 0, "Hs finite-dimensional, and 
Hs G Hs{'Hs) with maximal eigenvalue sz^ax o,nd mini- 
mal eigenvalue sz^in- For m > z^ax — Zmin 



V+{U)[is <E) P>„^] 



(B16) 



= [Is ® P>™-._+2_]F+([/)[l5 ® P>rn]- 

Hence, if cr (z S{'He) *s such that P>mO'P>m = cr, then 



ypeSiHs), yueUiHs). 



(BIT) 



Proof. Since m > Zmax — ^min it follows by equation (B13) 
that 



V+{U)[ls <S) P>m] 
N 
Yl X! \^n}{lp,i\Q\i'n'){lpn'\<»\l- Zn}{l~ Zn,\P>^. 
n.7i' — l 

(BI8) 



;> 



Let |a;) be an energy eigenstate of the reservoir. For 
(a;|y+([/)[ls ig) P>m] 7^ 0, a necessary condition is that 
X = I — Zn and / — z„' > m for some I, n, n' . We can thus 
write x-t-Zmax > x + Zn = I > m-\- Zn' > ?Ti-|-z,„i„. Hcncc, 

{x\V+{U)[ls «> P>,n] = for all X < m~ Z,nax + Zmin- 

It follows that [is ^ {Ie - P>m-z„,.^^+z„,J]V+{U)[is <S) 
P>m] = 0. This proves the lemma. D 
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The following Proposition tells us that a reservoir with 
a state that is L times removed from the ground state 
with respect to the maximal possible energy change D, 
can be used L times before there is any noticeable differ- 
ence in the set of channels that can be reached by using 
the reservoir. We refer to this as a quasi-catalytic prop- 
erty. Here XmaxiH) and Aniin(_ff) denote the maximal 
and minimal eigenvalues of the Hermitian operator H. 

Proposition 4 (Quasi-Catalytic states). Let S > 0, 
D (z N, and let Hsi i ■ ■ ■ , T~Lsl ^^ finite-dimensional, with 
Hsi G Hs{Hsi), such that K„ax{Hsi) - Xmin{Hs,) < sD 
for I = 1, . . . ,L. Suppose ct'-^' G S{He) is such that 



P>T,na^^^P~ 



>LD 



rW 



(B19) 



Define the sequence a^'"'^^' = A^ ^ (cr'-'-'), for I = 
1, ... ,i — 1, for pi e S{Hs,) and Ui £ U{Hs,). Then 

Cl""'^{a^^))=cl""'^{a^''^), / = l,...,i. (B20) 

Proof. By a repeated use of Lemma [8j it follows that 
P>iL-i+i)De^^'''P>iL-i+i)D = cr(''. Thus, by Lemma [t] 
we know that A+ ^^ (o-(')) = Ap,,t/, (o-*^''). By Lemmalrlwe 
can also conclude that C^''^^'' (cr(')) ^ C^i^Hs, (^rW). Fur- 
thermore, by Proposition[2]C'^''^s!(Ap,_^^[/,_j(cr('~i))) = 
C^'-^s, (■^(i-i)~) gy combining the above observations, 
we see that 

=C^-^-^<(A+_^,^,_^(a('-i))) 

The above reasoning can now be iterated, which yields 
the statement of the proposition. D 



3. Regenerative cycles 



Here we construct a protocol to maintain the set of 
reachable channels C+. Basically, what we do is to keep 
the state of the energy reservoir sufficiently far away from 
the ground state by pumping energy into it. 

The following lemma tells us that we can implement 

J. m 

the rigid translation a i— > A^crA^ within our model, 
if we have access to an ancillary system in a pure ex- 



cited state. By the observation that P. 



(Oi 



<)\J) 







for Zn 



> I > z„ 



and all j > we can directly use 



(B5| and (B6| to prove the following Lemma. 



Lemma 9. Let Ha G Hs^Ha) have the largest eigen- 
value sZjjjax with corresponding eigenstate lamax); o,i^d 
smallest eigenvalue s^min with corresponding eigenstate 



that 

Then 

A' 



Assume z^ax ^ z„ 

^ l^niax/ 



Let U £ UCHa) be such 

flmin). (B22) 



|a,„,„>(Qmax|.C/A 



(a) = A^ 



'aAl 



(B23) 



for all a G S{H]^) where A+ is as defined in equation 

(Bm. 



One may note that in this implementation there is no 
need to avoid any 'border zone'. 

In essence, the following proposition tells us that if the 
state of the energy reservoir has a support that is suffi- 
ciently far from the ground state, then we can maintain 
its power to induce channels on S (as characterized by 
the set Cj^ ®) merely by feeding energy into it between 
each use. 

Proposition 5 (Protocol for a catalytic cycle). Let 
s > 0, 13 G N, T-Lsi , 'Hs2 be finite- dimensional, and 
Hs, G H,{Hs,), Hs, G H,{Hs,), where X^^^AHs^) ~ 
Ainin(-ffsi) < sD and Xn^s.AHs2) - XminiHs^) < sD , with 
D G N. Let (Tin e S{'H%) be such that 



P 



>2D(7i-[iP>2D — Ci, 



(B24) 



For an arbitrary p\ G SCHsi ) o.nd an arbitrary Ui G 
UiHsJ, let 



Ki.uM'^)^ 



(B25) 



where the channel A'^ ^ is as defined in equation {BIO) 



Let T-La be two-dimensional, with Ha € Hs{'Ha) such 
that Ha has the largest eigenvalue sD with correspond- 
ing eigenvector |ai), and lowest eigenvalue with corre- 
sponding eigenvector |ao). Let Ua = |ai)(ao| -|- |ao)(ai|, 
and define 



A+ 

|ai)(ai|,;7A 



{^). 



(B26) 



where Af >, .jj is as defined in equation (BIO). Then 
" Kai„)=C+"''^~^(ai„) 



C 



— qS2,Hs 



(B27) 



P>2DContP>2D — CTout ■ 

Proof. Combining (B24) and Lemma Is] yields 
P>D^P>D = ^■ 



(B28) 



(B29) 



By further combining this with Lemma [o] gives ( B28| ) . 
Due to the latter we can, by Lemma [7] conclude that 



^S2,Hs2 f \ 

C, ((Tout) 



-'S2MS. 



(o^out)- 



(B30) 
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By (B29) and Lemma Is] it follows that (Tout = 
Ki}{ai\.UA^'^^ = ^\ai){a,\,UAi'^)- Hence Proposition g 
tells us that 



-'S2,Hs. 



(A|ai)(ai|,L'A(^)) =C 



S2,Hs^(- 



(^) 



(B31) 



Due to (B24) and Lemma^it follows that A^ ^ (cin) — 
Api,(7i(crin). By Proposition [2] 

C^-^-^(Ap„c;,(ain)) =C^-^-^(ai„). (B32) 



Furthermore, due to (|B24| it follows by Lemma [t] that 



C 



S2,Hs 



((Tin) =C_^_^' ''^(o-in) 



(B33) 



By combining the above chain of equalities, we obtain 
the proposition. D 



Appendix C: More general models 

In the previous sections we demonstrated that there 
in principle exist systems where coherence can be made 
catalytic. Here we briefly touch upon the question to 
what extent these phenomena survive in more general 
systems. 



to each other as ^0 increases (for a fixed L). One way 
in which this could happen if 1/'" would converge to a 
specific limit matrix as I — > +00. However, this is not 
the only possibility. For example, let H he a,n N x N 
Hermitian matrix, and define V^^'^ := exp{—if{l)H), 
where / : N ^ N is such that f{l + L) - f{l) -^ 
as I — )■ +00. [An example of this is f{x) — y^, 
for which \/x + L — ^x ~ Ll(2^fx)\ In this case 

(y('«), . . . , F('«+^)) ^i„_oo (^('°), . . . , 1^('°)). 

Hence, for a limited range of levels, the transformations 
become asymptotically uniform (in spite the fact that 
that they do not converge to a limiting operator). In 
Section |C2| we show that the Jaynes-Cummings model 
have this asymptotic uniformity. 

As a final remark we note that in Section |B] we could 
maintain the capacity of the coherence by injecting en- 
ergy into the reservoir. Furthermore, we did this solely 
'within' the model, i.e., we used the same class of uni- 
taries V^iJJ^ for the pumping, as for the operations on 
system S. It is a relevant question for which classes of 
these more general unitary operations that it is possi- 
ble (at least approximately). Unfortunately we do not 
provide an answer to this question here. 



2. The Jaynes-Cummings model 



1. General energy preserving unitary operations 

Here we consider Hamiltonians of the half-infinite lad- 
der model, but we relax the restriction on the dynamics, 
and only demand that it preserves the total energy. 

One can realize that a unitary operator V on 
the system and half-infinite ladder is block-diagonal 
with respect to the energy eigenspaces, i.e., V — 



p^)VP?\ if and only if 



N 



1^= E E KS'IV'n)(V'n'l®l^-^«)('-^«' 

^max-1 



(01) 



where each V^^^ = [% „']^„'=i is a unitary matrix. 

Recall that the harmonic oscillator model in Section iBl 
has the feature that it, sufficiently far from the ground 
state, essentially implements copies of a unitary interac- 
tion. We thus regain the model in Section |B] if we for 
all I choose F'^'' = y for a fixed unitary matrix V. This 
'uniformity' over the energy ladder appears to be closely 
connected to the catalytic properties of these models. It 
thus seems reasonable to suspect that in models that have 
some type of approximate uniformity, one should be able 
to find some approximate version of the type of behaviour 
that we have seen for the previous models. 

Suppose that the set of unitary matrices 
V^^°\ . . . ,V'^^°'^^^ would to become more 'similar' 



Here we numerically analyze the Jaynes-Cummings 
model [3 12] of a two-level system in resonant interac- 
tion with an electromagnetic field mode. As we shall see, 
numerical tests suggest that the higher the energy of a 
(suitably chosen) initial state of the reservoir, the longer 
its coherence survives sequential interactions with a col- 
lection of two-level systems. 

To this end, we consider the Jaynes-Cummings model 
[HH] for the interaction between a two-level system Hs — 
hi\i'o)i'o\ + h2\i'i)ipi\ and an electromagnetic field mode 
He = huja^'^a, with the interaction Hamiltonian 



^JC — g<^+ ® a + gcF- ® a) , 



(C2) 



where a, a) are the standard annihilation and creation 
operators [a^a^] = 1^;, a = X^z^i v ^1' ~ l)('l ^^"^ '^+ = 
|Vi)(V'o|,a_ = |V'o)(V'i|. 

We furthermore assume resonant conditions (i.e., the 
'matching' of the energy levels), so that hi — ho = hu, and 
choose hi = hujzi, ho = hujzo, with zo = 0,2:1 = 1. By 
standard textbook calculations, where we introduce the 
unit-free evolution parameter x := tg/H for the evolution 
time i, the evolution induced by Hjc can be written 



^-r.H,a^J2 E Vn,l\^n){i'n'\®\l-n){l-n'\ 
l—l n,n'— 

+ |V'o)(V^o|®|0)(0|, 

cos(xv^) — isin(xv^) 
-ism{x^/l) cos(xa/I) 



y« := 



15 



where we can write 



1 

1 



(C3) 



Hence, the dynamics in the JC-model provides an exam- 
ple of the type of asymptoticaUy uniform unitary opera- 
tors discussed in the previous section. The evolution on 
the reservoir is given by 



A.,|.0„)(^„|(a) :=Trs(e-"«|Vo)(V'o| ® ^e^"^) 

OO 

A^ :=^cos(a;V?) 10(^1 



B^ :=^sin(xv^)|Z-l)(Z|. 



1=1 



The sequential preparation procedure results in a evolu- 
tion process on the reservoir, which corresponds to an 
iteration of the channel A 



3:;\i'o)(M>o\ 



The nature of this 
dynamics is depicted in Figure |4] 

Here we attempt to mimic the sequential preparation 
procedure that we investigated in Section |A6| Again 
we thus use the (square of) the fidelity to measure how 
well the superposition |0) := (|V'o) — *|V'i))/'s/2 can be 
prepared from the ground state |'0o): by using the state 
a on the reservoir as a resource. More precisely, for a 
given X and a, our measure of 'quality' is 



,(a):=(0|$,,,(|^o)(^|)|0) 

^ OO 

: - + ^ sin(a:yr+I) cos(a;v^)Re(;|cr|; + 1) 



(C4) 



1=0 



Analogously to what we did in Section |A6| (where we 
optimized over all U) we should in principle optimize 
Qx{cr) over all x > 0. However, the dynamics of the 
JC-model can unfortunately be rather irregular (see e.g. 
[2 I40j ) why such an optimization may not be feasible. 
Instead we opt for a simpler (but possibly suboptimal) 
procedure, where we keep the value of x fixed (corre- 
sponding to one fixed unitary). For the sake of simplicity 
we choose x = 7r/4. For an initial state cr*^"^ we construct 



the sequence a^^~^^' = A, 



„|(cr^'^^) and the corre- 



7r/4,|i/>o>(Vi 

spending fidelities QT^/ii'^^''^)- For a suitably chosen Iq 

the initial state is u^^^ :— \'r]L,i(,){VL,io\- ^^ ^^ have seen, 
the evolution is asymptotically uniform. Hence, if L is 
small compared to ^Oj one can make the approximation 



0./4(a(°)) 



i+2(l-^)sin(vrVI/2). (C5) 



and thus obtain the (approximate) maxima 1 — 1/(2L) 
at 



lo{m) := {Am + l)"^ , m = l,2. 



(C6) 



Each graph in Figure [3] in the main text is obtain from 
an initial state that is a uniform superpositions stretching 
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FIG. 4: The state of the reservoir after iterated use. 

For an initial state cr'"' = \i]){r]\, with a very wide superpo- 
sition |r?) = E;^£o"i 0/^1000, the figure depicts \{l\a^^°°'>\l')\ 
(black corresponds to high values, and white low) where ct'^""' 



is the 100th iteration of a^''^'^^ = A, 



■/l.lV'o 



XMi^ 



W\ 



for 



X = 7r/4. The dotted lines corresponds to the values of I 
(and I') for which B^/i\l) = {B^/S') = 0). Since B^ in- 
duces loss of quanta from the reservoir, this suggests that 
weight will tend to accumulate at the crossings of the dotted 
lines. This figure shows an intermediate step in this process. 
The squares show the points lo{m) = (4m -I- 1)'^ around which 
(3^/4 takes a high value (for relatively narrow distributions). 
On these points we center the family of initial states that lead 
to the curves in Figure Isl (As opposed to the a^'^' that yields 
the present figure, those initial states only contains a narrow 
distribution of number states.) 



over 50 number states, centered around lo{m). Hence, 



JO) 



|^(m)^^^(m)|^ where |7?(")) 



Ef=-24l^o(™) 



A = - 
-n(m) 



Each graphs corresponds to the value F'^'"'' := Q.;r/4(o'm ) 
on the y-axis versus the number of iterations k on the x- 
axis. The different graphs correspond to the different 
values of m, where in the rightmost part of the figure, 
the curves correspond to m = 4, 6, 8, . . . , 28, 30 counted 
from below. These graphs suggest that the 'lifetime' of 
the coherence can be made longer, merely by increasing 
the average energy of the state. 

Since these graphs in Figure |3] are the result of an ap- 
proximate optimization, one should maybe not get too 
surprised by the fact that one see an initial increase of 
the fidelity in Figure [3| As seen, the largest increase hap- 
pens for the initial states with a relatively low number 
of quanta, which is where we indeed would expect the 
approximation to be the worst. The dotted line in Fig- 
ure |3] corresponds to the fidelity we would have obtained 
for the doubly-infinite ladder- model, i.e., 1 — 1/(2L) (see 
Section |A6| . Note that there is no reason to expect this 
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to be a bound for what can be achieved in the JC model. 
However, the gradual approach of the curves to this value 
suggests that the applied approximation is reasonable. 



One may wonder why we here use \'riL,io) &s ini- 
tial states, rather than coherent states, |a) := 
g-|a|V2 J2l^ jiJY72\l) IMSS], which would be the both 
convenient and conventional choice for this type of model. 
For the coherent states we could create a sequence of in- 
creasing average energy by increasing |a|. However, as we 
increase |a| we also increase the width of the superposi- 
tion over the energy eigenstates. Hence, in some intuitive 
sense the degree of coherence increases as we increase | a | , 
why it appears risky to use these states for the questions 
we investigate here. One should keep in mind though, 
that since we have not developed any operational mea- 
sure of the degree of coherence, discussions of this kind is 
a bit shaky. (For the doubly-infinite and the half-infinite 
ladder models, we could side-step the issue of measures 
of coherence by focusing on the set of induced channels 
Cia).) 



Appendix D: Prelude. Expected work-extraction 
without an energy reservoir 



Before we begin with the actual purpose of analyzing 
the role of coherence in expected work extraction by using 
explicit energy reservoirs (which we do in Section^ we 
first consider some simpler types of models that do not 
contain an explicit energy reservoir, and where one can 
show that the optimal expected work extraction (suitably 
defined) is given by the relative von Neumann entropy 
kTD{p\\GiHs)), where D{p\\a) = Tr(plnp) - Tr(plncr). 
The purpose of this exposition is partially to highlight 
the difference between using, and not using, an explicit 
energy reservoir, but also to gain some understanding 
for the ideas behind the (somewhat technically messy) 
derivations in Section |EJ On a technical level, the only 
things that we need from this section are Lemma [TT| and 
Corollary [2] 



In Section |D 1| we begin with a short 'prelude to 
the prelude', where we briefly recall the common ap- 
proach to implement expected work extraction using 
time-dependent Hamiltonians. More precisely, the sys- 
tem traverses through a path of Hamiltonians. These 
paths of Hamiltonians will implicitly play a role in the 
subsequent derivations, as we, in some sense, will simu- 
late such paths. In Section [D 2| we recall the concept of 
passivity, and the related setting for work extraction. In 
Section [D 3| we implement one version of the above men- 
tioned 'simulation', and show how this leads to optimal 
expected work extraction. 



Expected work extraction by varying the 
Hamiltonian 



A common approach (that can take many different 
guises on the level of assumptions and technical details) 
is to define expected work extraction in terms of changes 
of the Hamiltonian of the system. (For a handful of ex- 
amples, see e.g. ^01 1551 HTW5] . were [IS] provides an in- 
troduction.) In this setting we are given quantum states 
p and Hamiltonian H, and are allowed changes of Hamil- 
tonians H ^^ H' . The expected work gain of such a step 
is defined as 



W, 



yield 



{Hs ^ H's, p) := Tt{Hsp) - Tt{H'sp). (Dl) 



The contact with the heat bath is modeled as a replace- 
ment map p I— ?► ^h{p) '■— G{H). (As we will come 
back to below, there are variations on this, which can 
implement 'softer' types of thermalizations.) A process 
is defined as a sequence of shifts of Hamiltonians, sand- 
wiched by thermalizations, where the total work gain 
is defined as the sum of all the individual gains as in 
(Dl). This model is, for expected work extraction, a 



quantum generalization of the discrete classical model 
used in [Hj . (Unfortunately, this generalization does not 
make much sense as a generalization for e-deterministic 
work extraction.) Very much analogous to |14| . one can 
show that for a cyclic change of Hamiltonians, the op- 
timal expected work extraction in this setting becomes 
kTD{p\\G{Hs))- To see this, one can use the relation 
H ^ F{H) - kT\nG{H), where F{H) = -kT\nZ{H), 
and Z{H) ~ Tre^*^^, to show that the total work gain 
of an L-step cyclic process (TJ^^^ — _ff'^' — Hs), with 



P 



(0) 



p and p' 



C+i) 



$ 



_H-(i + 



1) (p*-'-*) becomes 



L-l 



WVieid = ^[Tr(iJ«pO - Tr(i/('+i)pO] 



1=0 



= ^D{p\\G{Hs)) - ^D{pL-i\\G{Hs)) 



L-l 



i^[l?(p,_i||G(i7W))-i?(pH|G(ffW))]. 

(D2) 



1=1 



Using <^h{G{H)) = G{H), together with the fact that 
the von-Neumann relative entropy is contractive un- 
der channels, i.e., D{r{p)\T{(T)) < D{p\a), it follows 
that D{pi\\G(H<-'))) = i?($H(„(p/-i)||$ff(o(G(iJ«)) < 
i^(p,_i||G(ff«)), and thus I^yieid < kTD{p\\G{Hs)) - 
kTD{pL-i\\G{Hs)). Due to the fact that D(-||-) > 0, it 
follows that within this model, for any cyclic process, we 
get the general bound Wyidd < kTD{p\\G{Hs)). 

As a side-remark one may note that, for the above 
derivation, we do not need to use the replacement map 
p !->■ G{H). As seen, we could, instead of ^h, use any 
channel that has G{H) as a fix-point, i.e., ^{G{H)) = 
G{H). (This implies a considerable freedom in imple- 
menting the softer types of thermalizations mentioned 
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above. However, we shall not need this freedom in 
our derivations, but will stay put with the, rather bru- 
tal, complete thermalization implemented by the replace- 
ment map.) 

The next question is if it is possible to saturate the 
bound VFyicid < kTD{p\\G{Hs)) within this model. A 
(limit) process that does achieve the bound can be de- 
scribed as follows: Change the initial Hamiltonian Hs 
into H', where the latter is such that p — G{H'). 
(Strictly speaking, this only works if p has full support. 
This can be treated via yet another limiting process, with 
a sequence of Gibbs states that converges to the state, 
akin to what was done in [M].) Next, we thermalize the 
system and perform an alternating sequence of thermal- 
izations and small changes of the Hamiltonians, keep- 
ing the eigenbasis fixed, until we reach a Hamiltonian 
H" that is iso-spectral to Hs- Finally, we perform an 
alternating sequence of thermalizations and small uni- 
tary transformations, until H" has been rotated back 
to Hs- The first of these steps has the expected work 

yield W^ieM = F{Hs) - F{H') + kTD{p\\G{Hs))- The 
second step, does in the limit of an infinitely fine dis- 
cretization, yield Wj^fjid = F{H') - F{H")- (Since we 
only have states that are diagonal in a fixed energy 
eigenbasis, the results for the ITR in [T3] are applica- 
ble.) For the final step of the protocol, let [/ be a uni- 
tary that transforms an eigenbasis of H" into an eigen- 
basis of Hs- We diagonalize U = X^n '2"^"lCn)(Cn|, 
for < 0,1 < 27r, and define the Hermitian operator 
A ■-— X]n<^nlC")(S"l- We can thus define a sequence of 
Hamiltonians if^^ := e-^^^^^H"e''^^/\ for / = 0, . . . , L, 
such that i?(°) = H" and i/*^' = Hs- Consider a pro- 
cess the sequentially changes i/'-'' to iJ*^'+^\ sandwiched 
by thermalizations. The work yield of this process is. 

wS^,^j:mHi^'>-H(^^^y)G{H(^'^)] 

=LTt[{H" - e-''^/^H"e''^/^)G{H")]- 

By inserting the expansion e~*'^/^if"e'^/^ = H" — 
i[A,H"]/L + 0{l/L?) into the above equation, one sees 
that W^vioid — J' as L — > oo. Hence, the total work yield 



yield 

W^yicld = <ield+ Al- 
dose to kTD{p\\G{Hs))- 



^yicid ean be made arbitrarily 



2. Work extraction and passive states 

The concept of passive states were introduced in the 
context of axiomatic characterizations of Gibbs states 
PSl Wf\ - See also recent results in [48j and ^49] which 
study the role of entanglement in this type of setting. 

We are given a quantum system S, with a Hamilto- 
nian Hs, and a quantum state p- We are allowed to 
perform arbitrary unitary operations (no matter whether 
they mix energy levels or not). The work gain of a uni- 



tary operation U is defined as 

Wyieid(C/, P, Hs) ■-= Tt{Hsp) - Tr{HsUpU^)- (D4) 

In light of the results in the previous sections, one may 
suspect that this construction implicitly assumes the ac- 
cess to ideal coherence. To analyze this issue in a quan- 
titative manner is precisely the purpose of Section [Ej 

A state p is passive with respect to a given Hamilto- 
nian Hs, if we cannot extract any energy from it in the 
sense of (|D4[ ), i.e., Wyieid{U,p,Hs) < for all unitary 
U [ini 147). it turns out [tSl gT] that a state is passive 
if and only if the state p and the Hamiltonian Hs has a 
common eigenbasis, and that in this basis the eigenval- 
ues of the density operator are ordered non-increasingly 
with increasing energy eigenvalues (or equivalently if the 
eigenvalues of p and G{Hs) co-decrease). 

It is straightforward to calculate the optimal extraction 
(e.g., by use of Lemma 12) 



roptiinal 
yield 



W^TJi™:= sup Wy^MU) 
ueuCHs) 



^Tt{HsPs) - inf TriHsUpsU^) 

ueuiUs) 

= -^D{ps\\G{Hs)) - -^D{X^{ps)\\XKG{Hs)))- 

(D5) 



■ \\j (p) means the eigenvalues of the op- 



Here At(p) > ••■-A, 

erator ps ordered non-increasingly. As seen from ( D5 ) 



we obtain the optimal work extraction precisely when we 
put the system in a passive state. One might wonder 
from where the /3 in ( D5 ) comes from, as we have not yet 
introduced any heat bath. At this stage it is merely a 
mathematical identity. [All the /3 on the right hand side 
of (D5) cancel, for any non-zero /3.] 



Due to the fact that relative entropy is non-negative 
^(•||') > 0, the optimal value can never be larger than 
kTD{p\\G{Hs))- If we somehow could 'arrange' for the 
term kTD{X^{p)\\X^{G{Hs))) to become zero, we would 
obtain kTD{p\\G{Hs)) as the optimum. 

An important observation is that a combination of 
passive states is not necessarily passive UHl SZI- As 
an example, for a passive state p, the tensor product 
p ® p may not necessarily be passive. For this we as- 
sume that the two subsystems have identical and non- 
interaction Hamiltonians, resulting in the total Hamilto- 
nian Hs-^ ®^s-2 + Isi ® Hs.2 ■ We furthermore allow arbi- 
trary unitary operations on the joint system "Hsi ® T~Ls2 ■ 
The Gibbs state G{H) have the special property that it 
is completely passive [461 I47j . meaning that it remains 
passive for arbitrary tensor products G{H)'^'^- 

Combining copies of a system is not the only way in 
which wc can change a passive state into a non-passive 
state. Another method is to add on a ancillary system 
yl in a Gibbs state G{Ha) for some Hamiltonian Ha- 
We furthermore assume that the ancillary system is non- 
interacting with S, such that the total Hamiltoni an is 
HsA ■■= i?s ® Ia + is ® Ha- Substituting 'S" in ([D5| 
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with SA with the new Hamiltonian Hsa, and the new 
initial state p(8) G'(i?^) [where we let the /3 in (D5 1 be the 



same as the (3 in the Gibbs state G{Ha) of the ancillary 
system] yields 

W^yTcL =Tr((i?5 + Ha)p ® G{Ha)) 

inf Tr((i/s + Ha)Up ® G{Ha)U^) 

UaUiHs'S'HA) 

^-^D{p\\G{Hs)) 

--^D{\\p®G{Ha))\\\KG{Hs)®G{Ha))). 

(D6) 

Comparing ( |D6[ ) with ( |D5| ) we see that the first term 
is unchanged [due to fact that D{p® G{Ha)\\G{Hs) ® 
G{Ha)) — D{p\\G{Hs))] while the second term is al- 
tered. lngeneviiiD{X^{p®G{HA))\\\^(G{Hs)®G[HA))) 
can be smaller than D{\^{p)\\\^{G{Hs)y), reflecting the 
fact that we can make a passive state non-passive by 
adding an ancillary Gibbs state. The purpose of the fol- 
lowing section is to show that the former term can be 
made arbitrarily small by choosing the ancillary Hamil- 
tonian Ha appropriately. 

The justification for the use of these ancillary Gibbs 
states is that we in principle always can obtain them by 
equilibrating the desired ancillary system with respect to 
the given heat bath of temperature T, i.e., an ancillary 
Gibbs state is a 'free resource' [9l[T0l[15]. More precisely, 
one can imagine that the ancilla system is put in contact 
with the heat bath. Next, we let it continuously equili- 
brate with the heat bath while we slowly de-connect it. 
(Since we here deal with questions about ideal thermo- 
dynamics, we do not worry about equilibration times.) 



3. Simulation of time-dependent Hamiltonians 
using a sequence of subsystems 

As is well known, relative entropy does not satisfy the 
triangle inequality. For the construction of our proofs, we 
make use of this 'failure' of the relative entropy. Lemma 
|11| shows that if we would regard the relative entropy as 
a (strange) distance measure, then the total sum of the 
'lengths' of the pieces of a path can be made arbitrar- 
ily small merely by making the division finer. (This is 
very closely related to the 'isothermal reversible paths' 
as described in [2].) 

Analogously to how we defined G{H), Z{H), F{H), 
etc., for operators, we here analogously define them for 
real vectors h = {hi,. . . ,hfq) e M^, by Gn[h) := 



-PK 



/Z{H), Z{h) := Y.r, 



'fShr. 



, F{h):=- j:\nZ (h) 



Lemma 10. Let h : [0, 1] — >■ M.^ have component-wise 
continuous second derivatives, then 



dx'^ 



F{h{x)) 



<max|/i"(a;)|-f 2/3max|/i',(a;)|^ (D7) 



By combining the above lemma with the Taylor ex- 
pansions, f{x) = f{xo) + {x — xo)f'{xQ) + r(x), with the 

error term in the Lagrange form r{x) :— ^^^^°' f"{^) for 
^ G [min(x, xo), max(a;, xo)] one can obtain the following 
lemma. 

Lemma 11. Let N E N and /3 > be given, and let 
h : [0, 1] — >■ M.-^ have component-wise continuous second 
derivatives. For each integer K > 2 

A'-l 



0<Y.D{G{h{^))\\G{hi'^)) 



< 



fc=0 



'kI^I 



K 



max max 

e[o,i] " 



\h'L{x)\ -\- B max max|/i'(a;)| 



(D8) 



In the following we let P( A^) denote the set of probabil- 
ity distributions over N symbols, and P+(A^) the subset 
of probability distributions with full support. 

Lemma 12 (Corollary IL4.4 on p. 49 in [SD]). Let a,b G 

R^^, then 



M 



M 



J2 "™^™ ^ Zl "' 



7n— 1 



M 

b < V a''' 6^ 

m—1 



(D9) 



Any permutation tt on the set {!,..., M} induces an 
operation H on P(M) by n^(p)m '■— P7r(m)' i-^-' ^^ P^"^" 
mute the elements in the index of p. A direct consequence 



of Lemma 12 is the following 



Corollary 2. Let q G P(A/) and p £ 

any permutation on {1. . . . , M}, then 

D{q^p^) < D{Il^iq)\\p). 



-{M), and vr be 



(DIO) 



Suppose that M := N^ , and suppose moreover that 
we choose to represent the elements of P(Af^) as L- 
dimensional tensors, i.e., as (pni,...,ni,)(ni,...,ni), where 
(ni,...,ni) G {1, . . . , A^}^^. We can obtain a special 
permutation on the index set {1, . . . , A^}^^ by permuting 
the 'index of the indices', i.e., (ni, . . . , n^) is mapped to 
(n7r(i), . . . ,«7r(L))j where it is a permutation of {1, . . . ,L}. 
This is of course merely a special case of a general permu- 
tation of {1, ... , A^}^. This permutation of the index of 
indices induces a corresponding operation 11 on P(A^^), 

by ll((,Pni,...,ni,jni,.. .,«£,; ~ (Pn„(i) ,...,«„(£,) jru ,...,nt ■ ^^ 

this is a special case of the more general permutation, 
it follows that Corollary [2] is applicable, which is impor- 
tant for the the proof of the following proposition. 

Lemma 13. Let N and q,p £ P+(A^) be given. For each 
integer K > 2 there exists a h G M such that 

D[{qG{h^)y\\{pG{h'')y) < Imaxlln^l^ (Dll) 

Proof. Since, q,p & IP-i-(A^) , it follows that log^n and 
log Pn are well defined and finite. Thus we can define 
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^n := — 4lng„, hl^ := — 4lnp„. We connect these two 
points with the path h : [0, 1] -^ R^ by h{x) := (1 



.K 



pN'' 



x)h^ + xh^ . Lemma 11 is apphcablc, and thus dehvers 
h''-'^ := h{^) with k = 0,l,...,K such that 

y DiGih'^-'') G(/i'=^^)') <-niax|ln^|2. (D12) 

^ \ ^ ' ' / K n p„ 



fc=0 

Define 

"■no,...,nK ■ "'no ' ^ ""n.K ' 



(no,...,nK) e{l,...,iV} 



K+l 



(D13) 



Hence /i^ e M^'^^', and 



G 



lO:K 



K:K\ 



no,ni ,...,n 



,(/.«) = G„„(/.'^-)...G„,(/i«-). (D14) 



On M we define the operation 11 by 

t,-^-^'^jns,no,ni,...,n/^_i ,n/^ • — ^no,ni ,...,n/^_i ,n/^ .n^ • \^^^j 

In other words, 11 is defined via a cychc permutation of 
the index of the indices of a. By Corollary [2] it follows 
that 

D[{qG{h^)4{pG{h^f 
<D(l{{qG{h^)) pG{h^) 
=D{Gih''-^)qG{h°-'') ■ ■ ■ G(/i^-i^^ 

pG{h^''^) ■ ■ ■ G{h^-^-'^)G{h'^-^) 
=D{q\\G{h'--^))+D{G{h^--^)\\p) 

+ J2D{G{h''-'')\\G{h''+'-'')). 



(D16) 



k=0 



By construction q = G{h ) and p — G{h 



,K:K\ 



(D12| the statement of the lemma follows. 



and by 

D 



Proposition 6. Let p e S+{Hs) and let Hs he a Her- 
mitian operator on lis- Then there exists a sequence of 
ancillary Hilhert spaces H^ and hermitian operators H^ 
such that 



w: 



opt 
yield 



lim 

X-f+oo 



Tr{{Hs + H^)p'E>G{H^)) 



inf TiiiHs + H^)Up G{H^)U^) 

UeU{'Hs»-H'^) 



^-D[p\\G{Hs)). 



(D17) 



Proof. Let q be the eigenvalues of p and p the eigen- 
values of G{Hs)- Since p G S+{'Hs) it follows that 
q G P+(A''). As p are the eigenvalues of a Gibbs state, 
it follows directly that p G P+(A^). By a repeated use 
of Lemma 13 we can construct a sequence (/i^)ifgN, 



where each h^ satisfies equation (Dll) in 



Lemma llSl Let {\a^)}^^i be an orthonormal b 

ar)(a; 



^ ^ ^ asis m 

an iV^^^<limensional Hilbert space H^ and define the 



can see that D{X^{p(E)G{H^))\\X^lG{Hs) (E)G{H^))) 



Hermitian operator H^ := X]fc=i 



^f 



One 



D{{qG{h^))^\\{pG{h^))^).' Co mbin ing this observa- 
tion with equations (D6| and ( |D11[ ), we obtain < 
kTD{p\\G{Hs)) - W% < 
proves the proposition. 



K 



|ln2iL|2 which 

n 



Appendix E: Expected work extraction w^ith an 
explicit energy reservoir 

In the context of passivity one defines (as we saw in 
Section D 2 ) the work gain as the decrease of internal 



energy of the system of interest under unitary opera- 
tions. Here we introduce an explicit energy reservoir, 
and thus model the system into which the extracted en- 
ergy is to be put. Consequently we define the expected 
work gain as the increase of expected energy in the reser- 
voir. Throughout this section we will exclusively make 
use of the model introduced in Section [X] 



1. Expected work extraction in fixed systems 

a. Operations on the system and energy reservoir 

Definition 2. Let s > 0, let Us be a finite- dimensional 
Hilbert space, with Hs G Hs{Hs), and let a G S{He)- 
Let szk and {ij^k) be eigenvalues and corresponding or- 
thonormal eigenvectors of Hs ■ Define 

N 

R^HQ):= Y. IV'n)(V'n|Q|^„')(V'n'|Tr(A^"-^"'a), 

n,n' — 1 

VQ G LiUs). 

(El) 

The channel i?^-^ should not be confused with ^^ jj 
introduced in Section fA 21 

Lemma 14. For s > 0, Hs G Hs{Us),and a G S^Ue) 

i?^ni) = i- (E2) 

In other words, the channel is unital and can thus not 
decrease the von Neumann entropy. 

So far in this investigations we have only dealt with 
expectation values of bounded operators, and have thus 
not needed to worry about the existence of these expecta- 
tion values. However, in this section we need to consider 
expressions like Ty{H^ a). Since H^ is an unbounded 
operator, we cannot guarantee that the expectation value 
is well defined for all elements a G S{'He)- One way of 
dealing with this would be to restrict S{T-Le) to elements 



20 



for which the product Hj^cr is a trace class operator {Q Let p G S{Hs) and a £ S*{He)- Then 



is trace class if Ti y/Q^Q < +00). In essence, we need 
a restriction such that {n\a\n') converges sufficiently fast 
to zero as n, n' — > ±00. However, for the purposes of this 
investigation there appears to be no strong reason to con- 
sider this technical issue in detail, why we in the follow- 
ing merely let S*{'He) denote some restriction of 5(H_e), 
where we (in one way or another) have made sure that all 
the relevant expectation values make sense. Note that all 
the operations we perform on the energy reservoir merely 
involve the interactions with a finite-dimensional space, 
why there is no reason to expect that these operations 
brings the state out of (a reasonably defined) S*{He)- 

Lemma 15. Let s > 0. Let T-Lg be a finite- dimensional 
Hilhert space, Hs G Hs{%s), P G S{'Hs), and a e 
S*{'He)- Then 



Tr 



{\is®H^E^]V{U)p®aV{U)'^ 



As) 



= Tr{H)^>a)+Tv{Hsp) 

-Ty(u^HsUR^^p)), 



(E3) 



for all U e ILJfHs) and all p G SlT-Ls): where V is defined 
as in Lemma\A and R^^ is as in Definition^ 



Proof. The 



first 

(s)i 

E 



part of the proof is to rewrite 

Tr( [iJs ® ifi + 
ls®H'^']p®a)-TT:(\Hs®lE]V{U)p®aV{Uy),^heTe 



Tr([l5 ® H]^']ViU)p (^ aViUY) 



we have used the fact that V{U) is energy conserving. 
By using equation ( |A5[ ) in Lemma [2| and use the fact 
that Hs is diagonal in the energy eigenbasis {|f/'n)}^=i, 
one can see that Tt{[Hs iE]ViU)p (g) aV{Uy) = 
Efc,,,pSZr(V'dC/V;')(^/'|C^|V'fe)(V'fc|p|V';)Tr(A^'=-^'a). 

The latter can be rewritten as Ti{W HsUR^'' (p)). D 



Given an Hermitian operator Q, we let X^{Q) denote 
the collection of eigenvalues of Q (counted with multi- 
plicities) ordered non-increasingly, while X^{Q) denotes 
the non-decreasing ordering. The following lemma is a 
direct consequence of Theorem 4.3.53 on p. 255 in jlSlj . 



sup TT{[is<»H'j^'']V{U)p<E)<jV{U)^) 
ueu{Hs) 

=Tr(i/i^V) 



+ iy^(p||Giii5Jj (E5) 

1 

1 

~~P 

Furthermore, there exists an element in UijHs) that 
achieves the above maximum, i.e., the 'sup' in the above 
equation can be replaced by 'max '. 



Dip\\G{Hs)) 
S{R^^p))-S{p) 
d(x^{R^Hp)) X^{G{Hs)) 



Compa red to (D5 



equation (E5). First 



we recognize two new features in 
ly, the term ^(i?^ ^ (p)) - S{p). Sec- 
ondly, that we have D{X^{R^^{p))\\X^{G{Hs))), rather 
than D{X^{p)\\G{Hs)). Both of these changes reflect the 
fact that the present m ode l explicitly takes into account 
coherence. Due to Eq. (E2), in Lemma 14 it follows that 



S{R^^{p)) > S{p). Furthermore by the properties of 
relative entropy £'(-||-) > 0. Hence, the work yield is 
bounded from above by kTD{p\\G{Hs))- As we shall see 
in Section 



E2 



the term S{R^^ (p)) — S{p) is determined 
by the relation between the degree of coherence in the en- 
ergy reservoir, and the degree to which the state p con- 
tains superposition between energy eigenspaces. Anal- 
ogously as to what we did in Section |D3[ the term 
D{X^{R^^ {p))\\X^{G{Hs))) can be made smaU by intro- 
ducing an ancillary Gibbs state. 



Proof. By combining Lemmas ^5] and 16 we obtain 



supc;Tr([l5 (g) H''^']VlJJ)p ® aV(Uy) = TriH^^'cr) + 
TriHsp) ~ Ek>liHs)Xi{R^Hp)). Next we use Hg = 
F{Hs) - kT\nG{Hs), and Xl{Hs) = F[Hs) ~ 
fcTlnAj^(G(i?5)). The rest of the proof is merely a rear- 
rangement to the terms. 

Finally we should show that the 'sup' can be replaced 
by 'max'. Let G\^{Hs) be the eigenvalues of G{Hs) or- 
dered non-increasingly, and let |^^) be corresponding 
eigenvectors. Let |0^J be eigenvectors of R^^{p) corre- 
sponding to the eigenvalues of X^{R^^{p)). The unitary 
operator U :— T^^^^^-i I'lpi) {(f'nl achieves the right hand 
side of equation 1E5|) . D 



Lemma 16. Let A, B be Hermitian operators on some 
finite- dimensional Hilbert space %. Then 



inf TtOJ^AUB) 

ueu(H) 



E^t(^)Al(i?)- 



(E4) 



Lemma 17. Let s,/3 > and M G N be given. Let Hs 
be a finite- dimensional Hilbert space, and Hs G Hs(Hs)- 



b. Adding a heat bath 

Analogously to what we did in Section |D 2[ we here 
model the heat bath by appending an ancillary system 
in a Gibbs state. 



As in Section D2 we assume that the system S and the 
ancilla A are non- interacting, i.e., the total Hamiltonian 
is HsA — Hs (g 1a + Is f^i Ha. (For the sake of notational 
simplicity we will in the following merely write Hs + Ha-) 
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We are allowed to freely choose the Hamiltonian Ha of 
the ancillary system, up to the condition that the eigen- 
values have to be multiples of the energy-level spacing s 
in the energy reservoir, i.e.. Ha G Hs{Ha)- We further- 
more assume that the ancillary system always starts in 
the Gibbs state G{Ha), and that the total initial state 
of the joint system is p G{Ha)- 

Lemma 18. 

Let s > and let Us o.'i^d T~La be finite- dimensional 
Hilbert spaces and Hg G HgiTis), Ha G HgiT-LA)- Let 
a e S*{'He), then 



R. 



Hs+Ha 



(p <E> G{Ha)) = R^Hp) <E> G{Ha), 



(E6) 



for every p G S{Hs)- 

Proof. We let hf^ — sz^ and ji/'f ) be the eigenvalues and 
eigenvectors oi Hs, and similarly let h^^ ~ sz^ and |V'm) 
be eigenvalues and eigenvectors of Ha- Hence, h^ + h^ 
and |i/'^,'0m) := \'^Pn)\i'm) ^i"6 the eigenvalues and eigen- 
vectors of Hs+Ha- The channel R^^^^^ can be written 
as 

N N 



n,n' — 1 m.,m' — l 

X Tr(A""-^"'- 



^'a] 



(E7) 



for aU Q G LIUs^Ra)- For Q := p(^G{Ha) we obtain 
the right hand side of (E6) by using the fact that G{Ha) 
is diagonal in {\'ip:^)},n- □ 

By a direct combination of Lemma |17| and |18| we obtain 
the following. 

Proposition 7. Let s > 0, Rs o.i^d Ra be finite- 
dimensional, and Hs € Hs{Rs), Ha € Hs{Ra)- Let 
p G S{Rs) and a G S*{Re)- Then 



supTr([l5«) 
u 

= Tr(4^V) 



He^]V{U)p (E) G{Ha) ® aViuY) 



D{p\\G{Hs)) 
SiR^Hp))-Sip) 



\d(x^{R^'{p) ® G{Ha)) \^{G{Hs) ® G{Ha)) 



P 



(E8) 



where the supremum is taken over all elements U G 

U{Rs)- 

Furthermore, there exists an element in U{Rs) that 
achieves the above equality, i.e., 'sup' can be replaced by 
'max '. 



In comparison with ( D6 1 we see that we again can af- 
fect the term D(A^(-)||A^(-)) by adding an ancillary Gibbs 
state, while the other terms (including S {R^ '=^ {p)) — S (p)) 
remain unaffected. 



2. Coherence in expected work extraction 



In light of Section A 3 it perhaps comes as no surprise 
that the level of coherence in the energy reservoir affects 
our ability to extract energy from the system. Here we 
first consider some natural special cases, and then pro- 
ceed to bound S{R^^ {p)) — S{p) in terms of the coherence 
in the energy reservoir. 



Diagonal system states and 
states 



energy-reservoir 



Let Hs = J2i ^iPi be the eigenvalue decomposition of 
Hs where Pi is the projector onto the eigenspace cor- 
responding to eigenvalue hi, where we let hi denote the 
distinct eigenvalues of Hs, i.e., hi ^ hi' ii I ^ V (as 
opposed to hn which is the complete list of eigenvalues 
including repetitions) . We define the following operation 



[Q]ff :=^PiQP,, VQgL(H). 



(E9) 



In words [Q]h removes the off-diagonal blocks of Q with 
respect to the cigenspaces of H (i.e., [■]h is what some- 
times is referred to as a 'pinching'). One can rewrite the 
channel R^^ in terms of the eigenprojectors Pi as 



i?fnQ)=E^'^^''Tr(A(^ 



i-h,,)/s 



a). 



LV 



(ElO) 



By using this reformulation of R^^ one can prove the 
following lemma. 

Lemma 19. If Q E ^{Rs) 's such that Q — [Q]hs then 

R^HQ) = Q, y<jeS{RE). (Ell) 

// a is such that a = [a] (s) , then 

R^HQ) = [Q]hs, yQeCiRs). (E12) 

By comparing Lemma [19] with Proposition [7] we can 
see that if p already is diagonal with respect to an en- 
ergy eigenbasis, then R^^ (p) = p and thus the term 
S{R^^{p)) — S(p) drops out. Furthermore, we obtain 
the following lemma, which tells us that if the state of 
the energy reservoir is diagonal, then the expected en- 
ergy yield can not depend on the off-diagonal elements 
of the initial state p. This implies that if the reservoir is 
diagonal, the expected work extraction can only depend 
on [p]i/s, which confirms the finding in |19j . 
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Corollary 3. Let s > 0, Hs and Ha be finite- 
dimensional and Hs G Hs{'Hs), Ha S Hs{'Ha)- Let 
p £ S{'Hs)- If cr & S*{'He) is such that a — [c]„(5}, then 

supTr([is ® i^ (g) H^j^''^]V{U)p (g) G{Ha) <E) aV{U)^) 
u 



/3 



D{[p]hs\\G{Hs)) 



-^D(^\^{[p]Hs<i^G{HA))\\x^{G{Hs)^G{HA))), 

(E13) 

where the supremum is taken over all elements U G 
UiHs + HA). 



b. Bounds on S(i?f -^ (p)) - S{p) 

As seen from the previous sections, the 'quahty' of the 
work extraction is partially determined by how much the 
effective channel R^^ increases the entropy of the initial 
state p. Here we determine bounds on this quantity in 
terms of the coherence in the energy reservoir. 

Proposition 8. Let S > and M e N. Let dimHs = 

N , and let Hg G Hs{T-Ls) with eigenvectors \^k) and cor- 
responding eigenvalues dzk- Let a G S^He)- Let 



the 



F:=[l-Tr(A^/-^^'Vr.,,^i, 



sup ||p-i?f-(p)||i<||F|| 



(E14) 



(E15) 



Proof. By using (El) in Definition [2] 



N 

p-R^^p)= J2 \^k){4'k\pHk'){i^k'\Fk^k'. (E16) 

k,k' = l 



N. We can thus without loss of generality assume that 

Due to the latter inequality we can conclude [e.g., 
via (A14|] that Fk,k' = \zk — Zk'\/L. By combining 
this observation with Proposition Is] and Lemma [4] we 
can conclude that \\p - R^-''{p)\\i < \\F\\ < iV(zmax - 

The following bound for for the von Neumann entropy 
is proved in |52j 

Lemma 21 ([SI]). Let dimH = N and p,a e S{H), 
then 

\S{p)-Sia)\<^\\p-a\\^logiN-l) + E[^\\p-a\\, 
S(x) :— —xhix — (1 — x) ln(l — x). 



(E18) 



Taking into account the fact that ^ is monotonically in- 
creasing on the interval [0, 1/2] and thus S(||/9 — cr||i/2) < 

S(7V(Zinax-2;min)/(2i)) if ||p-Cr||l < N {Zn^e,^- Zmin) / L < 

1 we obtain the following bound by combining Lemmas 
[20]and[2Tl 

Proposition 9. Let s > and let Hs be a Hilbert space 
with dimHs = ^ , and let Hs G HgiT-Ls)- Let z^ax be 
the largest eigenvalue of Hs divided by s, and let z,„in 
be the smallest eigenvalue of Hs divided by s. Assume 
L > A^(zmax - ^^min), and define 



<^L 



Then 



\VLJa)iVL,lol \VL,Io) 







Vl 



O/E- 



(E19) 



^fN{z^ 



+ 



-<'- 



'-NlogN 



2L 



(E20) 



with ^ as in LemmaWn 



If we define the matrix p = [{'4^k\p\'4'k')]kk'=i i* follows 
by equation ( |A17| ) that \\p - R^s[p)\\i = \\p * F\\i < 
||p||i||F|| = ll^^-^ D 

Lemma 20. Let "Hs be finite- dimensional with N = 
dim('H5) > 2, and let Hs G HsiHs), with eigenvalues 
szm- Let \r]L,io) := ^ J2i=o 1^ + ^o> G He, then 



sup Up 
pe5(«s) 



C)(...o|(^)ll^^ 



iV(z„ 



where 



max^^j^ Zk and Zmin 



L 



■ N 

mm^^j^Zk- 



(E17) 



Proof. In the case L < z^ax — 

true. This follows since both p and R, ^ ,, ,(p) are 

density operators and thus lip— ^i„^ \/„ i 



, Eq. (|E17|) is trivially 
x,,o|(^) are 
(P)ll <2< 



3. Achieving the 'standard' optimal expected work 
extraction 

Here we show that we can regain the 'standard' result 
on expected work extraction, namely that the maximum 
yield is kTD{p\\G{H)). We will obtain this as a hmit 
case, namely of a large density of states in the energy 
reservoir, high degree of coherence in the energy reser- 
voir, and very large (in Hilbert space dimension) ancil- 
lary systems. 

To do this, we introduce a number J that controls 
the density of states in the reservoir, compared to the 
energy levels in S. More precisely, if the system Hamil- 
tonian Hs has eigenvalues that are multiples of 5, then 
we assume that the energy spacings in the energy lad- 
der of the reservoir is s :— 5/ J. To provide coher- 
ence, we assume our standard uniform superpositions 
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\Vl,Io) = "121=0 I' + ^o)/VL, and thus get a number L 
that in some sense controls the degree of coherence in 
the reservoir. 

To construct the above mentioned hmit, we consider 
a fixed Hamihonian Hs S Hs{Hs), but a sequence of 
energy reservoirs, one for each K E N (where K is just a 
'dummy index'). We assume that the Kth reservoir has 
the energy spacing s{K) := 5/ J{K) and the initial state 
\'nL{K),io) (where Iq does not matter), with J{K) and 
L{K) being functions of K. We furthermore construct 
a sequence of ancillary systems A^ with Hilbcrt spaces 
Ti,^ with dimension N^^^, where N = diraHs- These 
ancillary systems are Gibbs distributed with respect to 
suitably chosen Hamiltonians. 

It maybe should be noted that in these derivations 
we limit ourselves to states p that have full support 
(p G S+{Hs))- This is essentially due to somewhat too 
weak bounds, which involves the logarithm of the small- 
est eigenvalue of p. Reasonably it should be possible to 
overcome this, but we will not do this here. 

Corollary 4. Let N e N, S,13 > be given. Let Hs 
be a Hubert space with divciHs = N, let p £ S-^-(Hs), 
and Hs G Hs{Hs)- Let L,J : N ^ N be such that 
\imK^+^[K/JiK)] = and\imK^+^[JiK)/L{K)] = 
0. Then there exist a sequence of ancillary Hilbert spaces 
■Hf with dim-Hf = 7V-^+\, and ijf e Hs/j(K){'Hf) 
such that 

lim WK.JiK)MK) = \d{p\\G{Hs)), (E21) 

A— >oo p 



a. Shifting the path to the grid 



where 



Wk.j{k).l(K) := sup Tr(i/^/^('^^a) 

-Tr{Hj^ '(Jl(k)) 



(E22) 



(E23) 



and where V{U) is as defined in (A5) for the reservoir 



Hamiltonian H^ , and 

f^L(K) ■^\VLiK),lo){VL{K),lo\, 

1 ''^^-' (E24) 

for some fixed Iq d 1,. 



This is a corollary of Proposition 11 The following 
subsections are devoted to the proof of this Proposition. 

It is maybe worth to point out that although the spe- 
cific construction we use to prove this proposition does 
require both increasing sizes of the ancillary systems, as 
well as increasing density of states in the energy reservoir, 
we do not strictly speaking prove that these are necessary 
conditions (although the author suspects them to be nec- 
essary in the generic case). Note also that this increase 
in density of states is in line with what was found in |19j . 
Recent results [53j may help to determine whether these 
are also necessary conditions. 



In Appendix |D 3| we demonstrated that one can ob- 
tain a discretized path of probability distributions, such 
that the sum of all consecutive relative entropies along 
the path can be made arbitrarily small. These prob- 
ability distributions we constructed as Gibbs distribu- 
tions G{h!^) of underlying energy level configurations h^ . 
The problem is that these energy level configurations 
may not necessarily satisfy our energy matching crite- 
rion, i.e., that each h^ should be a multiple of the energy 
reservoir spacing s. The purpose of this Section is to 
show that we can shift the energy level configurations h^ 
into new configurations h^ that are multiples of s, but 
at the same time obtain a sum of consecutive relative 
entropies Ef=o^ ^('^(^'')||<^(^''^^)) t^at is sufficiently 
close to the original Y.k~o D{G{h'')\\G{h^+^)). In M^, 
regarded as the space of energy level configurations h, 
the configurations that respect the energy matching form 
a square lattice. We construct the new configuration h 
by hn = s[/i„/sj, where [-J is the floor-function, which 
rounds of real numbers to the nearest smaller integer. 
(There is no particular reason to choose the floor-function 
in particular. We could equally well use [•] , or some other 
method to shift the original point h to somewhere close 
in the lattice.) 

Lemma 22. Let N,K E N and f3,s > be given. Let 
(/i^'OtLo with h'' e M^. Then there exists {h'')k=o ™*^^ 
h''^ E sZ^ , such that 

max \h'^ - h^J < s, D{G{h'')\\G{h'^)) < s/3, 

n=l,...,N (225) 

DiGih^mih")) < s/3, 
for k = 0, . . . , K, and such that 

Y^ D{GCh'')\\GCh''+^)) - Y^ D{G{h^)\\G{h''+^)) 



k=0 



fc=0 



< 2f3s{K + l) 

+ I3(e'''^ - 1)K max max \h':+^ - h'll. 

k=0,...,K~ln=l,...,N 



(E26) 



Proof For each k let h'^ := s[^h^\ , for n — 1, . . . ,N. By 
construction, this means that h'^ E sZ^, and 



|ft!;-^'j<5, /3s>ln5S>0, 



Z{h'^) 



^sp > Gnih) > ^-s0_ 



(E27) 



Gn{h) 
One can moreover realize that 

\Gn{h)-Gn{h)\<Gn{h){e'^-l). 



(E28) 



By Eq. (|E27| it follows that s/3 > D{G{h)\\G{h)). By 
the analogous reasoning one flnds sj3 > D{G(h)\\G{h)). 
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Observe that 

K-l 



^ D{G{h^)\\G{h''+^)) 



fc=0 



K-l N 



/3 E E^-e^')^^'^' - ^") + ^^Zih'') lnZ(fc") 

(E29) 



fc=0 n=l 



and the analogous statement for h. By making use of the 
inequahties in (E27I we find 



A'-l 



< 



K-l 



k=t) 



k=0 k 

K-1 N 

/3EEG"(^')(^^'-'^") 



k=0 n=l 
K~l N 



YY.Gn{h'){hl+^^hl) 



k=0 n=l 



K-1 N 



2/3s 



(E30) 



J2 E GnmCht' - hi+' -hi + hi) 

fc=0 n=l 
K-1 N 

fc=0 11=1 

+ 2/3s 
<2;3s(ii: + l) 

/f-l w 

+/3EE|g«(^')-^»(^')||'^^'-'^" 

k=Q n=l 



[By Eq. ( p28l )] 
<2Ps{K + 1) + ;9(e^^ - 1) E E G«(/i')|/ir' - K 



K~l N 



k=0 n=l 



<2/3s(is: + l)+/3(e"''- l)i\:max max \h: 

k—0 n—l,...,N 



(E31) 

D 



b. Simulation of the path 

Here we remind that P(A^) denotes the set of proba- 
bihty distributions over N symbols, and P+(A^) denotes 
the subset of probability distributions with full support. 

Proposition 10. Let N e N, (S,/3 > 0, andq,p£ P+{N) 
be given. For each K, J G N there exists a h^ e jZ^ , 



such that 



'^K\\i 



D[{qG{h^)Y {pG{h^)) 

1 K 

< 4/3(5-+ 2/35 y 






1 \ 1 ^'M 

1) max In — 

ri=l,...,JVl p„ 



(E32) 



1 h 9"|2 

— max In — . 



Proof. Since q,p E P+(iV), it follows that logg„ and 
logpn are well defined and finite. Thus we can define 
/i* := — 4lnq, h^ := — ilnp. We furthermore define 

h : [0, 1] -^ M^ by h{x) := (1 - x)h^ + xh^ for x G [0, 1]. 
The function h{x) = (1 — x)h'^ + xh^ clearly has a con- 
tinuous second derivative on each component. 
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Lemma 
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such that 



J in Lemma 
thus delivers h!"-^ := hf" , with h!'-'^ £ JZ^, 



For this function, let h :— h{-^), s 



< D{G{hiO))\\G{h°-^)) 



< 



< D{G{h^-'^)\\G{h{l))) < 



S(3_ 

T' 

(5/3 



(E33) 



and 



K-l 



Y, D{G{h^-^) G{h 



fe=0 
K-l 



-k+l:K^ 



YD{GiK^))Gihi^4l)) 



<2/3(5- 



fc=0 

.K+l 



(E34) 



J 



+ {e- 



K 



1) max I In — I 

'«=!,. ...at' pJ 



A combination of (E34| and Lemma 11 yields 



K-l 

Q <Y D(G{V'-'')\G{h''+^-'') 

k=0 
K-\ 

Y^(G(h^''^)\G(h''^^-'^ 

fc=0 



< 



fc=0 
K-\ 



(E35) 



E^(g(m^))||g(m^)) 

fc=0 



<2/3,5^^^^ + (eT _ 1) max | In ^ 

J n=l,...,JV pn 

+ - max |ln^|2 

K n=\,...,N pn 



By combining this with the two inequalities in Eq. ( E33 ) 
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we thus obtain 

0<D(G{h{0)) Gih"^-^] 



K-l 



■K^ 



G h 



fe=0 

+ D(G{h^-'')\\G{h{l)) 

1 K 

<Ap5j + 2P5 — 

+ ( e -f — 1 ) max In — 

n=l,...,N p„ 

+ -77 max In — . 

K n=l,...,N pn 



k+l-.K-^ 



(E36) 



c. Assembling it all 

A channel T : L{H) -^ L{H) is unital (or mixing en- 
hancing) if r(l) = 1. For a positive operator Q, let 
Amin(Q) denote the minimal eigenvalue. 

Lemma 23. // a channel T is unital then Aniin(r(/9)) > 
Amin(/o) for all p e S{'H). 

Proof. By [54j, or Theorem 2 in ^S], we know that for 
unital channels the output is more mixed than the in- 
put, i.e., X{T{p)) -< \{p). If the underlying Hilbert 
space has dimension N, this means Yl^l^i^k^ip)) — 
Sfc=i ^k(p) ^^^ "■ ~ 1, . . . , iV. Since F is a channel and 
thus trace preserving, we can conclude that Ainin(F(p)) = 



Define 



n-l ^If 



1 - ElZi At.(F(p)) > 1 - Efcli Kip) - A,„i„(p) 



n 



h^ ■= h°-^ + ... + h^ 

"'no....,nK ■ no ' '^ rij 

{no,...,nK) e{l,...,N} 



(E37) 



(E38) 



Since h''-^ G f Z^ it follows that 



Furthermore 

Gno,...,nK\'^ ) = tr„(,(« ■ ) • • • G,i^_j (n. ■ )G„^(/i ■ ). 

(E39) 
On M^ we define the operation H by 



(Ha) 



nsno,ni,...,n^_i,n/f 



ii,...,nR:_i ,nii:,ns • l^-Cj^Uj 



In other words, 11 is defined via a cyclic permutation of 
the index of the indices of a. By Corollary [2] it follows 
that 



D 
<D 



,^U 



(^iqG{h^)^\\{pG{h 
(^(gG(/l^'))||pG•(;^^) 
D(n(gG(/j°^^) • • • G(/i^-i^^)G(/i^^^)) 

pG(/iO^^)---G(/i^-i^^)G(/i^^^^) 
7^(G(/i^^^)gG(/i°^^') • • • G(/i^~i^^) 

pGih"-^) ■ ■ ■ G{h^-^-^)G(h'^-^) 
D{q\\G{h'-'<)) 

K-l 

+ J2 D{GCh''--'^)\\GCh''+^--'^)) 

fe=0 

+ D{GCh^-^)\\p). 



To avoid confusion in relation to Corollary HI one 
should note that in both Proposition [TT] and [121 J ^-nd L 
are not regarded as functions of K. 

Note furthermore the change of notation from R^^ to 

RJ ' ' ^ . This serves to underline the fact that rJ ' ^ 
depends on J, which should be kept in mind when us- 
ing Proposition |11| For a fixed density operator a — 
J2jj' ^ij'b)(/|j ^^ increase of J effectively means a 'de- 
crease' of the coherence in a compared to Hs ■ To see this, 
consider a transition from \tpn} to lipn')- The correspond- 
ing change in energy, s(z„' — z„), has to be compensated 
for by a change of energy s(z„ — z„/) = sJ{x„ — Xn') in the 
reservoir. The larger the J, the larger the compensating 
jumps in the reservoir have to be (counted in numbers of 
eigenstates) . In other words, an increase in J effectively 
decreases the 'width' of a in relation to Hs, and in this 
sense means a decrease in coherence. 

Proposition 11. Let N e N, S, 13 > be given. Let Hs 
be a Hilbert space with dim'Hs = N, let p G S-^-lHs), 
and Hs G HgiT-Ls), and let a e S*{He)- Then, for every 
J, K G N there exists an ancillary Hilbert space H^ with 
dim-Hf = N^+^, and H^ e Hs/j{Hf) such that 



(E41) 



<-^D{p\\G{Hs)) - ^ [s{Ri/'^-"-{p)) - S{p) 

1 K 

<A5-+25j 

+ ^(e^ - l)|ln[A„i„(p)A„,i„(G(iJs))]| 
+ ;^|ln[A„,i„(p)A„,i„(G(i/s))]|', 



Wj,k{ct) 



(E42) 



where 



- sup 



Tx(H^^I'^~g) 



By inserting ( |E35 ) we obtain the statement of the propo- 
sition. D 



Tr(i/f")a), 



(E43) 



~G ■.=Txs^a'^{V)p®G{hI)®oV{V)\ 
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yS/J.Hs 



and where V(U) is as defined in {AS) and Ra 

in DefinitionvA for the reservoir Hamiltonian H^ . 

^ ' ' ^ is unital (Lemma 14 1 we can, by 



Proof. Since Ra 

Lemma [231 conclude that 



X.nin{R'J''"Hp))>Xmin{p)- 



(E44) 



Hence, since p by assumption is lull rank, it follows that 
RJ ■ ' ^ (p) is full rank. Let q denote the eigenvalues of 

^(p))- Let p denote the 
X{G{Hs)). Since p by 



^S/J,Hsi 



R7 ■"■''' (p), i.e., q := X{R 
eigenvalues of G{Hs), i.e., p 

assumption is full rank, it follows that RJ ' ^ (p) is full 
rank, and thus q £ ¥+{N). Furthermore p G P+(Af) since 
G{Hs) is a Gibbs distribution. Furthermore, 



In- 



< - lnA,„i„(i?^/^^^-np)) - In X^UG{Hs)) 



By combining (E46), (E47), and (E49) we obtain the 
proposition. D 

As was mentioned above, an increase of J effectively 
means a decrease of the coherence in the reservoir state 
a. To guarantee that S{rJ' ' ^ {p)) — S{p) is small, we 
must thus choose the state a as to compensate for a large 
J. In the case a — \riL,io){VL,io\j "^^ must thus make sure 
that L:^ J. 

Proposition 12. Let N eN, 6,(3 > be given. Let Hs 

be a Hilbert space with diniHs = N , let p S S+{'Hs), and 

Hs E Hs{Hs)- With h^ the eigenvalues of Hs, we let 

Xn '■— hn/S, and Xmin := min„ a;„ and x„ 

Let J,K,L €N be such that L > NJ{xmi 

there exists an ancillary Hilbert space H^ with dimH^' 

N^+^, and H^ G Hs/jiV.^) such that 



— max„ Xn ■ 

'^ininj? llien 
K _ 

A — 



<-lnA,ni„(p)-lnA^i„(G(i/s)), 



(E45) 



where the last inequality follows by ( E44 ) 



Since q,p E P+(A^), we can conclude that the condi- 
tions of Proposition [TO] are satisfied, and we thus know 
that for each K, J there exists a h^ € 4^^ such that 



D({qG{h'')y {pG{h'')) 



K\\i 



1 K 

<4/35- + 2135 



+ (e 
1 



M 



J 

l)|ln[A„ 
max I InfAmi 



in(p)A,„in(G(i/s))]| 

{p)X^UG{Hs))]\\ 



(E46) 



where we have used (E45). Let Ti.^ be a Hilbert space 
with dimH^ = N^''^. We can construct a Hermitian 
operator H^ on V.^ that has h^ as its eigenvalues. By 
construction it follows that H^ g Hs/j{'H^). 
One can furthermore check that 



D[x^{Ri^J'"^p) ® G{H^))\\X^{G{Hs) ® G(i/f )) 
^D({qG{h^)y\\{pG{h^)y 



(E47) 

By assumption Hs G HsiHs), while iff G Hs/j{n^) 
and the energy reservoir has the energy spacing s = 6/ J. 
Note that 



HsiUs) C Hs/j{Hs) 



(E48) 



The assumption Hs G HsCHs) implies Hs G Hs/jCHs) 
and thus allows us to apply Proposition [7] on the state 



p (Xi G{H^) (X) a, which thus yields 



= '^Dip\lGiHs))^'/ 
1 



D{p\\G{Hs))^- SiR'J'^^np)) - S{p) 



/3 



:D 



(A^(i?^/^'^np) ® GiH^)) \\x^{G{Hs) ® G{H^)) 

(E49) 



0<^D{p\\GiHs))-WKj,L 



< 



J Xy] 



A* log A* + — ^ — 



P \L 



- PL 2 

1 K 

+ A5-+25j 

+ 4(e^ - 1)1 HXnun{p)X^UG{Hs) 

+ ^\HXmUp)Xm,n{G{Hs))]\\ 



where 



Wk. 



J,L 



sup 



TriH'J'a) 



s/j 



^Tviny-'aL) 



~TTsAViU)p(^ GiH^) (g> aLV{U) 



(E50) 



(E51) 



(E52) 



and where V{U) is as defined in lAw for the reservoir 



Hamiltonian H^ , and 



1 ^"^ 
V^ 1=0 



(E53) 



for some fixed l^ E'L. 



By letting J and L be functions of the parameter K, we 
can use Proposition [12] to prove Corollary |4] For this one 
can use the fact that due to limx->+oo[^(Ar)/L(Ar)] — 0, 
the inequality L{K) > A^J(A')(a;max — 2;nii,j) will hold for 
all sufficiently large K. 



Proof of Proposition 12 The eigenvalues ft.„ of Hs can 
be expressed in terms of multiples of the energy spacing 
s = 5/ J of the energy reservoir, as /if = z^S/ J. It proves 
convenient to re-express /i„ in terms of multiples in 5 
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(which we can do, due to the assumption 7^5 € Hs{T-Ls))- 
Hence, h^ — x„(5 and Zn = a;„ J, and thus Zmax = a^max^ 



anu ^niin ■^inin^- 



By assumption L > NJ{xyn_a->^ — a^min), and thus L > 
N{z-aiax — z-cnin)- Hcncc, Proposition ^ is apphcable and 
yields 



< S{Ry^'-"^{p))-S{p) 



(E54) 



< 



^inax ^miii 



2L 



NloeN + E 



^^ ( l-^inax ^Ta\n) 



2L 



max '^min 



< ^ ^ N\ogN + ^[ — 



and define Wj^k,l 



If we let a := ctl in Proposition 11 

Wj,k {o'l), and combine this witn (E54), we thus obtain 
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